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Chapter  1 
Introduction 


In  this  report,  the  interaction  of  electromagnetic  fields  with  objects  composed  of 
chiral  media  is  analyzed  using  eigenfunction  and  method  of  moments  techniques. 
The  term  chiral  media  was  first  used  by  Jaggard,  Mickelson,  and  Papas  in  1979 
[1],  who  defined  chiral  media  as  consisting  of  macroscopic  chiral  objects  randomly 
embedded  in  a  dielectric.  The  word  chiral  describes  something  that  is  handed,  i.e.,  an 
object  whose  mirror  image  can  not  be  produced  solely  by  rotating  and  translating 
the  original  object.  An  example  of  chiral  media  is  shown  in  Figure  1.1(a).  For 
comparison,  an  artificial  dielectric  is  shown  in  Figure  1.1(b). 

1.1  What  is  Chiral  Media 

Chiral  media  are  a  class  of  artificial  dielectrics  that  exhibit  optical  activity,  a  property 
of  many  biological  and  chemical  substances  which  occurs  naturally  almost  exclusively 
at  optical  frequencies  [2].  Optical  activity  is  a  reciprocal  property  characterized  by 
differing  wave  numbers  for  right  and  left  circularly  polarized  waves.  One  effect  is  to 
cause  a  linearly  or  elliptically  polarized  wave  to  rotate  its  axis  of  polarization  as  it 
propagates  through  the  optically  active  medium;  this  is  known  as  optical  rotation, 
circular  birefringence,  or  circular  double  refraction  [3].  The  other  effect,  circular 
dichroism,  is  the  difference  in  the  attenuation  rates  for  the  right  and  left  circularly 
polarized  waves,  thereby  causing  a  linearly  polarized  wave  to  degenerate  into  an 


1 


% 


^  § 

~  6s  < 


o°o° 

o°o° 


(a) 


(b) 


Figure  1.1:  Examples  of  (a)  a  chiral  medium  composed  of  a  three-dimensional  array 
of  randomly  orientated  helices  and  (b)  an  artificial  dielectric  composed  of  a  three- 
dimensional  array  of  spheres. 

elliptically  polarized  wave.  The  Faraday  effect  is  similar  to  optical  activity,  except 
the  latter  is  reciprocal  and  can  occur  in  isotropic  as  well  as  anisotropic  media,  while 
the  former  is  nonreciprocal  and  occurs  only  in  anisotropic  media. 

To  put  chiral  media  in  perspective  with  respect  to  other  media,  consider  the  con¬ 
stitutive  relationships  for  the  major  subdivisions  of  media  shown  in  Figure  1.2.  The 
constitutive  relationships  for  bianisotropic  media,  the  most  general  media,  are  shown 
at  the  top  of  the  figure,  where  the  explicit  dependence  of  the  constitutive  parameters 
on  space,  time,  etc.,  has  not  been  shown  for  clarity  [4,  5].  Research  on  bianisotropic 
media  can  be  applied  to  chiral  media  since  chiral  media  is  a  scalar  version  of  bian¬ 
isotropic  media  [6] — [18].  The  first  set  of  equations  shown  for  bianisotropic  media 
can  be  written  as 


(1.1) 


where  cD  and  H  have  the  same  units  (A/m),  E  and  cB  have  the  same  units  (V/m), 
and  in  combination  they  each  form  six-dimensional  vectors  (tensors  of  rank  one). 
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cD=P-E+L*cB 

or 

D  =  e  •  E  +  £  •  H 

H  =  M  E  +  Q  cB 

B  =  C-  E  +/Z-  H 

ANISOTROPIC 

BIISOTROPIC 

D  =  c  E 

D  =  eE  +  £H 

B  =  jl-  H 

B  =  (E  +  /iH 

ISOTROPIC 

CHIRAL 

D  =  €  E 

D  =  e  E  -  j  £c  B 

B  =  fiH 

Figure  1.2:  Constitutive  relationships  for  the  major  subdivisions  of  media. 


Of  particular  importance  is  that  all  of  the  constitutive  parameters  have  the  same 
units  and  together  form  a  six-dimensional  tensor  of  rank  two.  That  the  constitutive 
parameters  P  and  Q  individually  are  three-dimensional  tensors  of  rank  two  is  well 
documented  in  the  literature  on  bianisotropic  media;  however,  that  the  constitutive 
parameters  L  and  M  individually  are  pseudotensors  of  rank  two  is  often  neglected. 
This  characteristic  is  important  because  pseudotensors  are  not  invariant  with  respect 
to  improper  coordinate  transformations.  A  coordinate  transformation  is  improper 
if  the  determinant  of  its  Jacobian  is  negative;  one  example  is  the  transformation 
to  create  the  image  of  an  object  [6,  sec.  2.6,  6.2],  [19,  pg.  813].  This  translates 
into  a  minus  sign  as  was  shown  in  1971  by  Kong  when  developing  image  theory  for 
bianisotropic  media  [11].  A  similar  result  is  developed  in  Section  2.5  for  chiral  media 
using  the  volume  equivalence  theorem  of  Section  2.4.  The  same  result  can  also  be 
obtained  by  using  the  fact  that  (r  is  a  pseudoscalar  (pseudotensor  of  rank  zero). 

1.2  Overview 

The  remainder  of  this  chapter  presents  historical  background  on  the  constitutive 
relationships  for  chiral  and  optically  active  media,  as  well  as  an  outline  of  the  work 
on  chiral  media.  Chapter  2  presents  background  electromagnetic  theory  for  chiral 
media.  The  geometries  studied  in  this  report  are  shown  in  Figure  1.3.  Figure  1.3(a) 
shows  the  cylindrical  geometry  for  which  Bohren  [20]  developed  an  eigenfunction  so¬ 
lution  for  scattering  by  a  normally  incident  plane  wave.  This  solution  was  used  as  the 
base  for  developing  the  solution  for  scattering  by  the  multilayer  circular  chiral  cylin¬ 
der  presented  in  Chapter  3  and  shown  in  Figure  1.3(b).  The  eigenfunction  solution 
by  Bohren  was  also  used  as  the  reference  solution  when  developing  the  pulse-basis 
point-matching  method  of  moments  solution  for  scattering  by  an  inhomogeneous 
object  composed  of  chiral  media  in  free-space  presented  in  Chapter  4  and  shown  in 
Figure  1.3(c).  Chapter  4  also  presents  a  pulse-basis  point-matching  method  of  mo¬ 
ments  solution  for  scattering  by  an  inhomogeneous  object  composed  of  chiral  media 
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in  the  presence  of  a  perfectly  conducing  half-plane  as  shown  in  Figure  1.3(dj.  iiiis 
solution  could  be  generalized  for  scattering  in  a  half-space  and  verified  using  the 
image  theory  for  chiral  media  developed  in  Chapter  2  and  shown  in  Figure  1.3(e). 
Chapter  5  presents  a  spectral- domain  Galerkin  method  of  moments  solution  for  a 
microstrip  transmission  line  on  a  chiral  substrate  as  shown  in  Figure  1.3(f). 

1.3  History  of  Constitutive  Relationships  for  Op¬ 
tical  Activity 

This  section  presents  a  brief  history  of  the  constitutive  relationships  for  optically 
active  media.  In  the  fields  of  biology,  chemistry,  and  physics,  optical  activity  has  been 
studied  since  early  in  the  nineteenth  century  [2].  However,  it  was  not  until  nearly  a 
century  later  that  researchers  began  to  develop  theories  for  natural  optical  activity 
in  terms  of  electromagnetic  wave  theory.  The  early  theories  were  based  on  classical 
mechanics  [21,  pg.  616-635];  later,  quantum  mechanical  theories  were  developed 
[21,  pg.  703-723].  In  these  theories,  optical  activity  is  the  result  of  a  second-order 
term  in  the  dispersion  formula.  This  causes  a  time-harmonic  electric  field  to  create 
not  only  an  electric  dipole  but  also  a  linked  magnetic  dipole.  The  reverse  is  true  for 
a  time-harmonic  magnetic  field.  In  1937,  Condon  [3,  22]  summarized  the  quantum 
mechanical  theories  for  natural  isotropic  optical  activity  and  presented  constitutive 
relationships  of  the  form 


®  =  «*+/«-»* 

(1.2a) 

B  =  fio'H  +  f  S  +  g  ^ 

(1.2b) 

where  g  is  the  rotatory  parameter  and  /  is  small  with  only  a  second  order  effect  on 
the  phase  velocity  and  no  effect  on  optical  rotation. 
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(e)  Image  Theory 


(f)  Microstrip 


Figure  1.3:  Overview  of  geometries  studied  i 
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this  report. 


In  1959,  Moiiitt  and  Moscowitz  extended  this  work  with  constitutive  relationships 
of  the  form 


D(w,t)  =  e{u)  E(u;,f)  +  /(u;)  H(u;,f)  -  g{u)  —  H(u;,f)  (1.3a) 

B(w,<)  =  g(u)H(u,t)  + f{u)E(u,t)-\- g(u>)  ^-E{u>,t)  (1.3b) 

where  the  time-dependence  was  not  dropped  [23,  pg.  652].  They  showed  that  the 
real  and  imaginary  components  of  the  constitutive  parameters  e(u>),  fi(w),  /(<*>),  and 
g( u>)  are  related  by  the  Kramers-Kronig  relations  [5,  6] 


c»  -  -§/ 


ds 


<(s) 

a 2  —  u)2 


ds 


(1.4a) 

(1.4b) 


where  ^(a>)  =  ('{u)  —  j("{v)  and  (  —  e,  fi,  /,  g.  Moffitt  and  Moscowitz  showed  that 
the  optical  rotation  is  proportional  to  g'  and  the  circular  dichroism  is  proportional  to 
g";  therefore,  these  two  measured  quantities  are  related  by  Kramers-Kronig  relations. 
This  is  true  because  /  is  small  and  affects  only  the  average  phase  velocity  with  a 
second  order  effect.  This  can  be  seen  by  examining  the  wave  numbers  for  right  (k/f) 
and  left  (k/J  circularly  polarized  waves  propagating  in  an  optically  active  media 
described  by  Equation  (1.3): 


kn 


kL 


2 

v  9 
c 
2 

u  9 
c 


(1.5a) 


(1.5b) 


The  factor  kn  —  k[ ,  contains  only  g ,  and  it  is  this  factor  that  is  responsible  for  optical 
rotation  and  circular  dichroism.  Additional  information  on  optical  activity  and  how 
the  constitutive  equations  are  derived  can  be  found  in  the  previously  mentioned 
references  and  [24] — [36]. 
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1.4  Microwave  Models  for  Optical  Activity 

The  optical  activity  of  helical  molecules  led  researchers  to  consider  large-scale  mod¬ 
els  at  microwave  frequencies.  Some  of  the  earliest  experimental  work  on  this  topic 
was  done  in  1920  and  1922  by  Lindman  [37,  38],  who  measured  the  rotation  of 
microwaves  propagating  through  randomly  oriented  copper  helices.  In  1956,  Win¬ 
kler  [39]  repeated  Lindman’s  work  and  obtained  slightly  different  results.  Similar 
measurements  were  done  in  1957  by  Tinoco  and  Freeman  [40]  for  oriented  arrays  of 
helices.  They  showed  that  the  array  of  helices  rotated  the  polarization  of  incident 
plane  waves  over  a  wide  frequency  range.  The  strongest  optical  rotation  and  cir¬ 
cular  dichroism  corresponded  to  the  case  of  A  =  2 L/k,  where  k  =  1,2,3,...  and 
L  is  the  length  of  the  wire  in  the  helix.  However,  these  researchers  were  interested 
in  how  their  results  matched  models  for  optical  activity,  not  in  the  electromagnetic 
properties  of  arrays  of  wire  helices. 

1.5  Electromagnetic  Scattering  and  Optical  Ac¬ 
tivity 

One  of  the  first  researchers  to  study  scattering  by  optically  active  bodies  using 
electromagnetic  wave  theory  was  Bohren;  a  previous  researcher  used  the  approximate 
Mie  scattering  technique  [41].  Bohren  used  constitutive  relationships  of  the  form 
[42,  43] 

D  =  eE  +  cteVxE  (1.6a) 

B  =  iiH  +  ^VxH  (1.6b) 

where  a  and  (3  are  the  parameters  which  account  for  the  optical  activity.  The  pa¬ 
rameters  q  and  (3  are  equal  in  the  absence  of  an  externally  applied  magnetic  field. 
His  work  on  scattering  by  electromagnetic  waves  was  limited  to  eigenfunction  solu¬ 
tions  for  homogeneous  spheres  [42],  spherical  shells  [44],  and  cylinders  [20].  He  also 
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presented  an  approximate  solution  lor  arbitrarily  shaped  inhomogeneous  L-alies  us¬ 
ing  the  Rayleigh-Debye  approximation,  which  is  valid  only  for  small  bodies  with  low 
density  [45].  However,  his  interest  in  optically  active  bodies  was  based  on  modeling 
the  measurement  of  biological  particles  such  as  viruses  [46]. 

1.6  History  of  Chiral  Media 

The  history  of  chiral  media  starts  in  1979  with  the  introductory  paper  by  Jaggard, 
Mickelson,  and  Papas  [1],  They  examined  the  approximate  scattering  by  a  single 
perfectly  conducting  short  wire  helix  and  concluded  that  a  randomly  oriented  col¬ 
lection  of  such  objects  would  behave  macroscopically  as  a  medium  with  constitutive 
relationships  of  the  form 

D  =  eE-y£cB  (1.7a) 

H  =  -B  -  j£cE  (1.7b) 

where  £,  is  the  chirality  admittance  and  the  e^1  time  convention  is  used.  The  next 
paper  on  chiral  media  was  not  published  until  1982  and  considered  the  transition 
radiation  generated  by  a  charged  particle  passing  through  a  plate  composed  of  chiral 
media  [47].  Later  research  has  extended  this  work  to  Cerenkov  radiation  in  infinite 
chiral  media  [48,  49]. 

Chiral  Interfaces  and  Slabs 

Silverman,  in  1985,  examined  scattering  by  an  achiral/chiral  interface,  where  the 
chiral  medium  was  described  by  either  the  Condon  (Equation  (1.2)  with  /  =  0)  or 
the  Born  (Equation  (1.6)  with  /3  =  0)  forms  of  the  constitutive  relationships  [50]. 
Additional  research  has  been  done  on  scattering  by  achiral/chiral  interfaces  [51] — 
[57],  chiral/chiral  interfaces  [58],  the  interface  of  mirror-conjugated  chiral  media  [59], 
and  periodic  achiral /chiral  interfaces  [60].  Researchers  have  also  studied  scattering 
by  chiral  slabs  [54,  61]  and  chiral  slabs  backed  by  perfect  conductors  [62] — [64]. 
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Related  research  has  studied  wave  propagation  along  a  chiral  slab  [65].  Creation  of 
an  anisotropic  chiral  slab  using  slabs  of  uniaxial  media  has  also  been  studied  [66]. 

Scattering  by  Chiral  Bodies 

Research  on  scattering  by  objects  composed  of  chiral  media  has  expanded  since  the 
early  work  by  Bohren  [20],  [42] — [45].  In  1985,  Lakhtakia,  V.  K.  Varadan,  and  V.  V. 
Varadan  presented  results  for  scattering  by  spheroidal  objects  using  the  T-matrix 
method  [67].  Later  work  by  Uslenghi  extended  the  spherical  problem  to  a  single 
chiral  layer  on  a  sphere  wTith  an  impedance  surface  [68].  Related  work  has  been  done 
on  the  eigenmodes  of  a  perfectly  conducting  sphere  filled  wdth  chiral  media  [69]. 
Research  by  the  author,  presented  in  Chapter  3,  extended  the  homogeneous  circular 
cylinder  studied  by  Bohren  [20]  to  a  multilayer  circular  chiral  cylinder  with  the  center 
cylinder  described  by  an  impedance  surface  [70].  Additional  work  by  the  author, 
presented  in  Chapter  4,  extended  scattering  to  objects  of  arbitrary  cross-section 
using  the  volume  equivalence  theorem  [71,  72].  Numerical  results  were  presented 
for  two-dimensional  objects  in  free  space  [71,  72]  and  in  the  presence  of  a  perfectly 
conducting  half-plane  [73].  An  alternative  approach  requiring  fewer  unknowns  has 
been  developed  by  Rojas  [74]. 

Basic  Research 

The  majority  of  the  early  research  concentrated  on  the  electromagnetic  properties 
of  infinite  chiral  media  [75] — [87]  and  the  Dyadic  Green’s  functions  for  one,  two, 
and  three  dimensions  [88] — [90].  There  is  the  possibility  that  this  research  may  be 
applied  to  modeling  vegetation  layers  for  remote  sensing  [88].  This  research  has  been 
applied  to  antennas  radiating  in  infinite  chiral  media  [91,  92]  and  in  the  presence  of  a 
chiral  sphere  [93].  Related  research  has  been  done  on  a  point  dipole  radiating  inside 
a  chiral  sphere  [94].  Surface  integral  equations  have  been  developed  for  scattering 
by  perfect  electrically  conducting  bodies  in  an  infinite  chiral  medium  [95]. 
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Construction  of  Chiral  Media 


Research  into  the  construction  of  microwave  chiral  media  has  been  limited  to  the¬ 
oretical  studies  of  scattering  by  beaded  helices  [96,  97]  and  measurements  of  chiral 
media  constructed  with  wire  helices  embedded  in  a  dielectric  [61,  98,  99].  Related 
research  has  studied  the  mixing  of  small  chiral  bodies  in  an  infinite  chiral  media 
[100,  101]. 

Applied  Research 

In  1988,  an  excellent  introductory  article  on  chiral  media  was  published  by  Engheta 
and  Jaggard  [102].  Since  then,  research  has  expanded  into  parallel  plate  waveguides 
[103],  circular  waveguides  [104] — [109],  periodic  structures  [110],  lenses  [111],  and 
microstrip  antennas  on  chiral  substrates  [112,  113]. 
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Chapter  2 
Theory 


This  chapter  presents  electromagnetics  theory  for  the  analysis  of  the  interaction  of 
electromagnetic  fields  with  objects  composed  of  chiral  media.  Section  2.1  presents 
the  constitutive  relationships  and  basic  notation  for  chiral  media.  Section  2.2 
presents  a  technique  for  transforming  the  coupled  wave  equations  of  chiral  media 
into  a  set  of  uncoupled  wave  equations  in  terms  of  right  and  left  circularly  polarized 
waves.  Section  2.3  presents  circular  vector  potentials  which  have  properties  similar 
to  those  of  the  usual  magnetic  and  electric  vector  potentials,  except  that  they  result 
in  circular  rather  than  linearly  polarized  fields.  This  property  makes  them  useful 
for  field  expansions  in  chiral  media.  Section  2.4  presents  the  volume  equivalence 
theorem  for  chiral  media,  w'hich  is  used  in  Chapter  4  for  a  method  of  moments  cal¬ 
culation  of  scattering  by  arbitrarily  shaped  bodies.  The  volume  equivalence  theorem 
is  also  used  in  Section  2.5  to  develop  image  theory  for  a  chiral  body  over  a  perfectly 
conducting  electric  or  magnetic  ground  plane. 

2.1  Constitutive  Relationships 

This  section  presents  the  constitutive  relationships  and  notation  used  in  this  report. 
All  fields  and  currents  are  considered  to  be  time  harmonic  with  the  e]ul1  time  depen¬ 
dence  suppressed.  The  constitutive  relationships  for  a  chiral  or  an  isotropic  optically 
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active  medium  [6,  sec.  8.0)  can  be  written  as 


(2.1«) 

(2.1b) 


D  =  eE  -  j(r B 
H  =  -B-jtrE 

fl 

where  pL  is  the  permeability,  e  is  the  permittivity,  and  the  pseudoscalar  £r  is  the 
chirality  admittance  of  the  medium.  The  medium  is  lossy  if  /r,  e,  or  £c  are  complex.  If 
£r  =  0,  then  Equation  (2.1)  reduces  to  the  constitutive  relationships  for  an  isotropic 
achiral  medium.  To  simplify  the  following  developments,  Equation  (2.1)  can  be 
written  as 


D  =  erE  -  (2.2a) 

B  =  fiVL  +  j/i&E  (2.2b) 

where  the  effective  permittivity,  ec,  of  the  chiral  medium  is  defined  as 

er  =  e  +  ix£c2 .  (2-3) 

2.2  Chiral  Wave  Equation 

This  section  presents  the  wave  equations  for  chiral  media  and  a  technique  for  trans¬ 
forming  the  coupled  E  and  H  wave  equations  into  uncoupled  wave  equations  in 
terms  of  right  and  left  circularly  po':  ized  waves,  which  can  be  solved  by  the  classi¬ 
cal  approach. 

Regardless  of  the  medium,  the  fields  (E,  B,  D,  H),  produced  by  the  impressed 
electric  and  magnetic  currents  (J1,  M'),  are  related  by  Maxwell’s  equations  [114], 


-  V  x  E 

=  ju  B  +  M' 

(2.4a) 

V  x  H 

=  ju>D  +  J' 

(2.4b) 

V  •  D 

“  P 

(2.4c) 

V  •  B 

=  0 

(2.4d) 
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where  p  is  the  electric  charge  density  given  by  p  —  j'V  •  J/u>. 

Now  consider  a  chiral  medium  with  constitutive  parameters  (//,  e,  (,  ).  The  fields 
(E,  H)  produced  by  (J1,  M')  in  the  chiral  medium  a^e  obtained  by  substituting 
Equation  (2.2)  into  Maxwell’s  equations  resulting  in 


-  V  x  E 

=  jupH  —  u;/i£rE  4-  M' 

(2.5a) 

V  x  H 

=  E  -f  H  J 

(2.5b) 

V  -E 

=  pl* 

(2.5c) 

VH 

=  ~HcP/€. 

(2.5d) 

Equation  (2.5)  can  be  written  in  matrix  form  as 


(2.6a) 


where 


1*1  = 


U)pt,c  -jup 

juier  u>p.(c 


Using  Equation  (2.6),  the  source-free  wave  equation  is 


(2.6L) 

(2.7) 

(2.8) 


Following  the  work  of  Bohren  [42,  43],  although  for  constitutive  relationships  differ¬ 
ent  than  those  used  here,  the  coupling  caused  by  \K]  in  the  wave  equation  can  be 
removed  by  diagonalizing  [A-]  such  that 

k/i  0 

\K  ]  =  \A\  W'.  (2.9) 

0  -kL 

Using  linear  algebra  techniques,  [A]  is  found  to  be 

[A]=  1  1  (2.10) 

j/V r  ~j  Nr 
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with  the  chiral  wave  impedance  given  by 


and  the  chiral  wave  numbers  given  by 


=  u}y/pTc  ±  a )fl(c. 


(2.11) 


(2.12) 


Substituting  Equation  (2.9)  into  Equation  (2,8)  and  multiplying  through  from 
the  left  by  [/t]-1  results  in 


U]-'V2 


To  simplify  Equation  (2.13)  define 


kn  0 

0  -kL 


(2.13) 


(  H  j  =  1  e"  j 

where  E/<  and  E/  can  be  shown  to  be  the  electric  fields  of  right  and  left  circularly  po¬ 
larized  waves  with  wave  numbers  kn  and  ki,  respectively  [43].  Then,  the  uncoupled 
source-free  wave  equation  in  chiral  media  is 


k2RER 


(2.15) 


Multiplying  Equation  (2.14)  through  from  the  left  by  [y4]  gives 

/  E  \  T  1  1  1  /  E«  \  /s 

=  •  (2-16) 

\  H  /  j/vc  ~j/Vc  \  El 

Substituting  Equation  (2.16)  into  Equation  (2.6)  results  in  the  uncoupled  Maxwell’s 
equations  for  chiral  media 


\  1  1  ji i 


1  -JVc 


\  ~kL El  ) 

P  1  ~  VrCr 
2e  1  +  Vc(c 


(2.17a) 


(2.17b) 
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From  Equation  (2.16),  the  magnetic  fields  are  given  by 


(2.18) 

The  total  electric  or  magnetic  field  is  the  sum  of  the  right  and  left  circularly  polarized 
electric  or  magnetic  fields 


E  =  Eft  +  E^ 
H  =  Hft  +  HL. 


(2.19a) 

(2.19b) 


2.3  Circular  Vector  Potentials 

This  section  develops  the  right  and  left  circular  vector  potentials,  which  simplify  the 
field  expansions  for  chiral  media.  These  potentials  are  analogous  to  the  magnetic 
(A)  and  electric  (F)  vector  potentials.  In  an  homogeneous  source-free  achiral  region 
the  magnetic  vector  potential  A  is  obtained  by  solving  the  wave  equation 

V2A  +  k2A  =  0  (2.20) 

where  V  x  A  =  B,  V  •  A  =  —jiv/xe^,  and  $  is  the  electric  scalar  potential  [115]. 
However,  in  an  homogeneous  source-free  chiral  region,  the  corresponding  wave  equa¬ 
tion  is 

V2A  +  2u;/z£cV  X  A  +  k2A  =  0  (2.21) 

using  the  same  definitions  as  used  for  Equation  (2.20)  and  the  constitutive  rela¬ 
tionships  for  chiral  media.  Solutions  for  equations  of  this  form  are  known  [78,  88]; 
however,  an  alternate  solution  Is  presented  below  based  on  the  same  method  that 
was  used  in  the  previous  section  to  uncouple  the  wave  equation  for  chiral  media. 

Following  the  techniques  used  for  vector  wave  functions  [114,  sec.  7.1]  the  right 
(R)  and  left  (L)  circular  vector  potentials  are  defined  as 

R  =  aip(kn)  (2.22a) 

L  =  aip{kL)  (2.22b) 
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where  a  is  an  arbitrary  constant  unit  vector  and  rp(k)  is  a  solution  of  the  scalar  wave 
equation 

v2x(j{k)  +  fcV(fc)  =  o.  (2.23) 

From  these  vector  potentials  the  right  and  left  circularly  polarized  electric  fields  are 
formed  using 

E R  =  V  X  (r+^-Vxr)  (2.24a) 

\  / 

EL  =  Vx  |l  -  y  V  x  l)  .  (2.24b) 

The  right  (or  left)  circular  vector  potential  component  Ry  (or  Lv )  produces  a  right 
(or  left)  circular  to  y  field  RCj  (or  LCj  ),  just  as  the  magnetic  vector  potential 
component  Ay  produces  a  transverse  magnetic  to  y  field  TMy. 

2.4  Chiral  Volume  Equivalence  Theorem 

This  section  develops  the  volume  equivalence  theorem  for  chiral  media  [114,  sec.  1.6], 
[115,  sec.  3-11],  [116],  [117],  [118,  p.  327].  The  volume  equivalence  theorem  allows 
a  chiral  scatterer  to  be  replaced  by  free  space  and  equivalent  electric  and  magnetic 
volume  polarization  currents  (J,  M).  From  Maxwell’s  equations  the  incident  or 
free  space  fields,  denoted  (E1,  H1),  of  the  impressed  electric  and  magnetic  currents 
(J1,  M' )  are  related  by 

—  V  X  E'  =  jufi0ll'  +  M'  (2.25a) 

V  x  H'  =  ju>e0 E*  +  J*  (2.25b) 

where  (/z0,  eo)  are  the  constitutive  parameters  of  free  space. 

Now  consider  an  inhomogeneous  chiral  medium  with  constitutive  parameters 
(/x,  e,  £r).  The  total  fields,  denoted  (E,  H),  produced  by  (J‘,  M')  in  the  inhomoge¬ 
neous  chiral  medium  are  given  by  Equation  (2.5)  and  can  be  written  as 

—V  x  E  =  jiMfioli  +  —  n o)H  —  u>n£cE  +  M'  (2.26a) 

V  X  H  =  jwe uE  +  ju)(ec  —  e0)E  -I-  u>fi(c H  -f  J'  (2.26b) 


where  juifiuH  and  ju>e0E  have  been  added  and  subtracted  from  the  right  hand  sides 
of  Equations  (2.26a)  and  (2.26b),  respectively,  to  simplify  the  next  step. 

By  definition  the  scattered  fields,  denoted  (Es,  Hs),  are  the  difference  between 
the  total  fields  and  the  incident  fields, 

E5  =  E  -  E'  (2.27a) 

Hs  =  H  -  H'.  (2. 27b) 

Subtracting  Equation  (2.25a)  from  (2.26a),  Equation  (2.25b)  from  (2.26b),  and  mak¬ 
ing  use  of  Equation  (2.27)  yields 

-  V  x  F5  =  jujfiuHs  -f  [ju(fi  -  /i0)H  -  ufi(cE]  ■  (2.28a) 

V  X  H5  =  jW uEs  -f  [ju)(ec  -  e0)E  -f  u>/z£.H]  (2.28b) 

By  comparing  Equation  (2.28a)  to  (2.25a)  and  Equation  (2.28b)  to  (2.25b),  it 
can  be  seen  that  the  scattered  fields  appear  to  be  produced  by  the  equivalent  electric 
and  magnetic  volume  polarization  currents 

J  =  ju>(ec  -  e0)E  +  u>/x£c H  (2.29a) 

M  =  ju/(/i  —  //u)H  —  ufi(cE  (2.29b) 

radiating  in  free  space.  Note  that  (J,  M)  are  nonzero  only  where  the  parameters  of 
the  inhomogeneous  chiral  medium  differ  from  those  of  free  space. 

2.5  Image  Theory  for  Chiral  Media 

This  section  develops  image  theory  for  a  chiral  body  over  an  infinite  perfectly  con¬ 
ducting  electric  or  magnetic  ground  plane  [119].  Image  theory  defines  an  equivalent 
problem  for  currents  and  material  bodies  over  an  infinite  ground  plane,  in  which  the 
ground  plane  is  removed  and  replaced  by  the  image  of  the  currents  and  material 
bodies. 
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A  clurai  medium  was  originally  defined  as  consisting  of  macroscopic  chiral  objects 
randomly  embedded  in  a  dielectric  [1].  From  this  definition,  the  image  of  a  material 
body  consisting  of  chiral  media  does  not  have  the  same  material  parameters  as  the 
original  body,  since  the  image  body  consists  of  macroscopic  chiral  objects  of  the 
opposi+e  handedness. 

The  chiral  volume  equivalence  theorem,  developed  in  Section  2.4,  is  used  to 
replace  the  chiral  body  by  free  space  and  equivalent  electric  and  magnetic  currents. 
Next,  conventional  image  theory  is  used  to  replace  the  ground  plane  by  conventional 
image  currents.  Finally,  the  chiral  volume  equivalence  theorem  is  used  to  replace 
the  image  currents  by  the  image  of  the  original  chiral  body. 

Figure  2.1(a)  shows  the  original  problem  consisting  of  the  impressed  currents 
producing  the  total  fields  (E,H)  in  the  presence  of  a  chiral  body  with 
material  parameters  (/i,e,£c)  over  an  infinite  perfectly  conducting  electric  ground 
plane  at  y  =  0.  For  simplicity,  the  ambient  medium  is  considered  to  be  free  space 
with  material  parameters  (/q,,e u).  In  Figure  2.1(b),  the  chiral  volume  equivalence 
theorem  is  used  to  replace  the  chiral  body  by  free  space  and  the  equivalent  electric 
and  magnetic  volume  polarization  currents  given  by  Equation  (2.29). 

In  Figure  2.1(d),  conventional  image  theory  for  currents  is  used  to  replace  the 
ground  plane  by  (J'^M"),  the  image  of  the  impressed  currents,  and  (J',M'),  the 
image  of  the  equivalent  currents,  which  are  related  to  the  total  fields  (E,H)  in  the 
upper  half  space  by 

J'n,v(xi-y,z)  =  Tiw(ec  _  €o)En,\(x,y,  z)  u>fi(cHH.\  {z,y,  z)  (2.30a) 

~y,z)  =  ±ju>(fi  -  (x,y,z)  ^  ufi£cEi-i'\  (x,y,z)  (2.30b) 

where  the  subscripts  H  and  V  refer  to  the  horizontal  (*,£:)  and  vertical  (y)  vector 
components,  respectively.  The  image  fields  below  y  =  0  are  given  by 

E'n,\\x^-y^z)  =  TEh.\ {x,y,z) 

HH,v(xi-yiz)  =  ±H,pV(x,y,z) 


19 


(2.31a) 

(2.31b) 
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Figure  2.1:  Image  theory  for  a  chiral  body  over  a  perfect  electric  conductor  de¬ 
rived  using  the  chiral  volume  equivalence  theorem  and  conventional  image  theory 
for  currents. 
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Figure  2.2:  Image  theory  for  a  chiral  body  over  a  perfect  magnetic  conductor. 

Combining  Equation  (2.30)  and  Equation  (2.31),  (J',M')  are  related  to  (E',H')  in 
the  lower  half  space  by 

J'(x,-y,z)  =  ju>{ec- e{))E'(x,-y,z)~u;nZcH'(x,-y,z)  (2.32a) 

=  ju>(n  -  nu)ll'(x,-y,z)  +  u>n(cE\x,  -y,z).  (2.32b) 

Finally,  comparing  Equation  (2.29)  with  Equation  (2.32)  shows  that  the  image  of 
the  chiral  body  has  material  parameters  of  (/x,e,  —  £c),  as  shown  in  Figure  2.1(c). 
Similarly,  as  illustrated  in  Figure  2.2,  the  image  of  a  chiral  body  over  a  perfectly 
conducting  magnetic  ground  plane  also  has  material  parameters  of  (/r,e,—  £c). 

This  section  has  developed  image  theory  for  a  chiral  body  over  an  infinite  per¬ 
fectly  conducting  electric  or  magnetic  ground  plane.  In  both  cases  it  was  shown  that 
the  image  of  a  chiral  body  has  the  same  material  parameters  as  the  original  body, 
except  that  the  chirality  admittance  of  the  image  is  the  negative  of  the  original.  In 
essence  the  mirror  image  of  a  chiral  object  is  that  object,  such  as  a  helix,  with  the 
opposite  handedness. 
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Alternatively,  the  same  results  can  be  obtained  using  tensor  analysis  and  the  fact 
that  £ r  is  a  pseudoscalar  (pseudotensor  of  rank  zero).  Pseudoscalars,  unlike  scalars 
(tensors  of  rank  zero),  are  not  invariant  with  respect  to  coordinate  transformations. 
The  image  of  the  chiral  body  is  obtained  by  spatial  inversion,  an  improper  coordinate 
transformation,  which  introduces  the  sign  change  [6,  sec.  2.6,  6.2,  8.3],  [19,  pg.  813]. 


Chapter  3 

Scattering  by  a  Multilayer 
Cylinder 


This  chapter  presents  an  efficient  recursive  eigenfunction  solution  for  scattering  by 
a  multilayer  circular  chiral  cylinder,  with  or  without  a  surface  impedance  center 
cylinder,  for  transverse  magnetic  (TM)  and  transverse  electric  (TE)  normal  incidence 
plane  waves,  as  shown  in  Figure  3.1.  The  special  cases  of  zero  or  infinite  surface 
impedance  correspond  to  a  perfect  electric  or  magnetic  conducting  center  cylinder, 
respectively. 

The  solution  for  scattering  by  an  M  layer  chiral  cylinder  requires  finding  an 
eigenfunction  expansion  in  each  layer,  then  matching  continuity  of  tangential  electric 
and  magnetic  fields  at  the  M  interfaces.  The  M  =  1  or  homogeneous  cylinder  has 
been  analyzed  by  Bohren  [20].  In  that  case,  solving  directly  for  the  four  unknown 
coefficients  is  manageable.  However,  if  even  one  more  layer  is  added  the  coefficients 
in  the  eigenfunction  expansion  are  algebraically  complex  and  unwieldy  to  find.  This 
can  be  seen  in  Appendix  A,  which  presents  an  eigenfunction  solution  for  scattering 
by  a  perfectly  conducting  cylinder  coated  by  a  single  homogeneous  layer  of  chiral 
media.  An  alternate  technique  is  to  set  up  and  numerically  solve  a  4A/  x  4 M 
matrix  equation  which  enforces  boundary  conditions  at  the  M  interfaces.  To  avoid 
the  computer  CPU  time  and  storage  associated  with  large  matrix  equations,  here 
we  generalize  Richmond’s  method  [120,  121]  to  apply  to  multilayer  circular  chiral 
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Figure  3.1:  Geometry 
impedance  center  cylinder. 
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cylinders,  The  present  technique  requires  t lie  multiplication  of  M  4  x  4  matrices 

to  determine  the  coefficients  in  the  eigenfunction  expansion.  The  difference  between 
the  solutions  for  TM  and  TE  incident  plane  waves  is  extremely  simple,  differing  by 
only  a  plus/minus  sign  at  one  step  in  the  solution.  Further,  the  addition  of  a  surface 
impedance  center  cylinder  involves  only  one  additional  matrix  multiplication. 

3.1  The  Eigenfunction  Expansion 

In  this  section,  the  eigenfunction  expansion  for  a  general  multilayer  circular  chi¬ 
ral  cylinder,  with  or  without  a  surface  impedance  center  cylinder,  is  presented  by 
expanding  the  electric  and  magnetic  fields  in  terms  of  right  and  left  circularly  polar¬ 
ized  waves  using  the  notation  presented  in  Section  2.2.  The  geometry  of  the  cylinder 
is  shown  in  Figure  3.1,  where  the  constitutive  parameters  and  the  outer  radius  of 
layer  m  are  (/zm,  em,  £m)  and  rm,  where  m  =  1,2,..., A/.  The  right  and  left  wave 
numbers  of  the  mth  layer,  fcj?  and  A:™,  and  the  wave  impedance,  7 are  defined  by 
Equations  (2.12)  and  (2.11).  Note  that  m  =  1  is  the  innermost  material  layer  and 
m  =  M  is  the  outermost  material  layer.  The  outer  radius  of  the  impenetrable  surface 
impedance  center  cylinder,  when  included,  is  denoted  tv  The  external  medium  is 
free  space,  with  parameters  (/v,  eu),  wave  number  ky  =  and  characteristic 

impedance  77,1  =  yj /qj / Cu ■ 

First,  consider  the  fields  external  to  the  cylinder.  These  fields  may  be  expanded 
as  an  infinite  sum  of  vector  wave  functions,  Mn  and  Nn,  which  are  related  by  [114, 
sec.  7.1]: 

VxN„  =  fcMn  (3.1a) 

V  xMn  =  &Nn.  (3.1b) 

These  vector  wave  functions  are  solutions  to  V  X  V  X  C  —  fc2C  =  0  and  V  •  C  =  0, 
where  C  =  Mn  or  Nn.  In  cylindrical  coordinates,  they  may  be  written  as 

NW(fc)  =  ze}n4,Z^\kp)  (3.2a) 


25 


(3.2b) 


M !r>(k)  =  j,^*Z^(k,)  -  4  e’"*Z^\kp) 

where  Z ^  is  a  Bessel  function  of  type  p  ( Z ^  =  Jn,  Z ^  —  Y„,  Z\^  =  H^\ 
Z^  =  H and  the  prime  indicates  the  derivative  with  respect  to  the  argument. 

In  the  case  of  a  TM^  plane  wave  normally  incident  from  <f>  =  180°,  the  incident 
fields  may  be  written  as  [115,  sec.  5-8] 


OC 


E'  -  ze~]k°x  =  £  r"N(’)(fco) 

n=  — o c 

(3.3a) 

H'  =  -y-e-’k°x  =  £  rnMi'\k0) 

Vo  Vo„  =  -oc 

(3.3b) 

where  Equation  (3.3b)  is  obtained  from  Equation  (3.3a)  using  H  =  xE  and 
V  X  Nn  =  ^Mn.  For  the  TE^  case,  to  maintain  the  symmetry  of  the  solution,  the 
incident  fields  are 


OC 


E-  =  yje~*  0X=  ^  J-"M<')(fco) 

n  =  —  oc 

(3.4a) 

H1  =  z—e~jk<}X  =  -  £  i-nN(’>(fc„). 

Vo  Vo„^-oc 

(3.4b) 

Due  to  the  chiral  cylinder  the  scattered  fields  will  have  TM^ 

and  TEz  compo- 

nents;  therefore,  the  scattered  fields  are  expanded  as 

ES  =  £  r"  [onNMM  +  fc.M| W(*b)] 

n  =  — oc 

(3.5a) 

Hs  =  -  £  j~n  kM[^)  +  6nN(-')(*o)] . 

n=  —  oc 

(3.5b) 

In  the  chiral  cylinder,  the  eigenfunction  expansion  must  account  for  the  rotation 
of  polarization  inherent  to  chiral  media.  In  a  z-independent  problem,  this  results 
in  a  coupling  between  the  TM?  and  TE^  fields,  which  prevents  the  eigenfunction 
expansion  from  being  written  as  a  simple  superposition  of  TMz  and  TE7  fields,  as 
was  done  in  Equation  (3.5)  for  the  scattered  fields  in  the  free  space  region.  How¬ 
ever,  the  eigenfunction  expansion  can  be  written  as  a  superposition  of  right  and  left 
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circularly  polarized  fields.  This  is  done  by  combining  the  vector  wave  function,  Mn 
and  Nn ,  to  form  right  and  left  circularly  polarized  vector  wave  functions  [20] 


e£  =  MW(*„)  +  N!r>(*„)  (3.6a) 

E  &  =  (3.6b) 

which  satisfy  the  vector  wave  equation  (2.15)  and  Maxwell’s  equations  (2.17a)  and 
(2.17b)  in  a  source-free  chiral  region.  Then,  the  fields  in  layer  m  of  the  chiral  cylinder 
may  be  represented  as 

Em  =  E  riCEl'U^Ei1’  +/:E'j,>„+s:Ei,>]  (3.7a) 

n=  — oc 

H”  =  —  E  j-  [CEi't  -  CEi’l  4  rC  -  jJEfi]  (3.7b) 

n=  —  oc 

where  the  c™,  d™,  f™ ,  and  g™  are  unknown  coefficients  and  Equation  (3.7b)  is 
obtained  from  Equation  (3.7a)  using  Equation  (2.18)  of  Section  2.2.  Note  that 
the  fields  in  each  layer  of  the  chiral  cylinder  are  expanded  in  terms  of  inward 
and  outward  E(  l)  propagating  right  E n  and  left  El  circularly  polarized  waves,  thus 
requiring  four  coefficients  per  layer. 

Expanding  Equation  (3.7)  to  get  an  explicit  representation  for  the  field  compo¬ 
nents  in  each  layer  produces: 

e."  =  E  i~"  [cM^rp)  -  CJnKp) 

Tl—  —  OC 

-  sTWOT/O]  e’"*  (3-8*) 

HT  =  —  E  r"  [c Mk%p)  4 

+CH!?\Kp)  +  £  »?>(*"<>)]  *•"*  (3.8b) 

e;  =  -  E  j-'  [c Wp)  +  <K(Kp) 

n=  — oc 

4 £HW(KRp)  4  gZH^Wp)]  <4"*  (3.8c) 

=  —  E  i-  «(*M  -  <K(K p) 

7]rn  n—  _  oc 
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+/r  *:'”(*»/>) -»r»n*z>)]  r 
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(3.8f) 


3.2  Determining  the  Coefficients 


This  section  presents  a  recursive  technique  for  determining  the  coefficients  of  the 
eigenfunction  expansions  presented  in  the  previous  section.  These  unknowns  in¬ 
clude  c"1,  d”' ,  /"’,  and  g”‘  in  each  of  the  M  layers  as  well  as  the  au  and  bn  from 
Equation  (3.5).  This  technique  is  a  generalization  of  that  developed  by  Richmond 
for  achiral  multilayer  circular  dielectric/ferrite  cylinders  [120,  121].  The  unknown 
coefficients  are  determined  by  applying  the  boundary  conditions  of  continuous  tan¬ 
gential  electric  and  magnetic  fields  at  the  boundary  of  each  layer,  which  can  be 
written  as 


£yn  + 1 

II 

JS1 

“  3 

(3.9a) 

pm  +  l 

~  L<t> 

(3.9b) 

jjrn  +  l 

II 

N  3 

(3.9c) 

rjm  +  1 
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1? 

1! 

(3.9d) 

Using  Equation  (3.8),  the  boundary  conditions  can  be  expressed  as 
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MWp) 

-UKp) 

WKKp) 

-h^Kp) 
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J'nWp) 

WV'RP) 

RL!),(Kp) 

^MkRP) 

±r\krp) 

tWnP) 
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±WWp) 

-±W(k?p) 

(3.10) 


(3-11) 


The  resulting  relationship  between  the  coefficients  for  layers  m  and  m  +  1  can  be 
represented  in  matrix  form  by 


rm  +  \ 

’  c ' 

dm+l 

n 

fm  + 1 

J  n 

=  M 

< 

fn 

.  sr] . 

.  3n  . 

where 


[ym]  =  [zm+1]  1  [ zm ] . 

After  simplification1  the  matrix  [Tm]  is  given  by 


vL'4\RY 


U<‘-4>(L) 
-X<‘-‘>  (L,R) 
-ui'-'W 


U^(R) 

X^(R,L ) 

-vl!'\R) 

-X<<’>(R,l) 


X  «4J,(L,fl) 
U^(L) 
-X^(L,R) 


(312) 


(3.13) 


(3.14) 


’The  symbolic  algebra  program  MACSYMA™  supplied  by  Symbolics,  Inc.  of  Cambridge, 
Mass,  greatly  aided  in  this  simplification. 
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where 


U^(S)  = 


x 'S™\S,T)  = 


yks+-rm  [—  +  l]  •  [zW(kZrm)ZM'{k2+'rm) 

^  7]m 

-ZP'(k?rm)ZM(kZ*'r„)} 

—  L,m  +  l  [ 

3- -  r-~  —  -1  '  [Z<p>(fc?rm)Zfo>'(fc?+Vm) 

4  T1m  L 


(3.15a) 


(3.15b) 


and  S ,  T  are  L  or  i?. 

3.2.1  Perfectly  Conducting  Center  Cylinder 

If  the  innermost  cylinder  is  a  perfect  electrical  conductor,  the  boundary  condition 
of  zero  tangential  electric  fields  at  p  =  r0  can  be  written  in  matrix  form  as 


n!?HK ,r„)  -irMH'-i  /,! 

ffi!),(*M  af'Wn.)  s; 


~J„(k'nr„)  J„(k'br„) 

-MKr u)  J  [  < 

using  Equations  (3.8a)  and  (3.8c). 

Solving  Equation  (3.16)  for  f]n  and  g\  in  terms  of  c^  and  d]n  produces 

[  n  1  \  y  °  1  f  c\  1 


(3.16) 


v} i 


(3.17) 


where 


I?,  =  +  (3.18a) 

lAn 

=  s&s;  ‘318b> 

131  -  (3-18c) 

=  -4-  [UK'«) W(k'Hv0)  +  J:(klr„)ff<’>(klr„)]  (3.18(1) 


(3.18b) 

(3.18c) 
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and 


4.  =  jr„)  +  Hm'(k<Rr„)Hm(k'LT0). 


(3.19) 


To  define  a  [Fu]  that  can  be  used  with  Equation  (3.12),  Equation  (3.17)  can  be 
written  as 


(3.20) 


V 

V 

d} 

=  M 

d' 

n 

n 

f'n 

0 

.  9n  . 

0 

where 


i 

0 

0 

0 

0 

1 

0 

0 

yo 

*31 

y° 

■*  32 

0 

0 

yo 

1  u 

y° 

M2 

0 

0 

r  0  10  0 

[y°  I  =  .  (3.2i) 

*■  J  I/O  t  m  r\ 


3.2.2  Surface  Impedance  Center  Cylinder 

If  the  innermost  cylinder  is  described  by  TM  and  TE  surface  impedances  the  bound¬ 
ary  conditions  at  p  =  r0  are  given  by  [122] 


z  _  E* 
ZTE  ~  - H J 

Z  -  E‘ 

Ztm  -  T, 


(3.22a) 

(3.22b) 


where  the  perfect  electric  conductor  in  Section  3.2.1  is  the  special  case  of  Zje  —  0 
and  Zjm  =  0. 

Using  Equation  (3.8)  these  boundary  conditions  can  be  expressed  as 


Zr£-  +  &HW(k'RT „)  +  gW^klrJ  = 

V\  1  J 

k  +  i'J’M  r„)  +  11  ffW(k’R  r„)  + 

Ztm^  [c’-WJr’v)  -  i'Mkl r„)  +  fW{k'HT0)  -  jW'Wu)]  = 


(3.23a) 


[c'„Uk'Rr„)  -  KMk',r„)  +  fW’Hkk r„)  -  fli JVi1' (*i.r0)]  .  (3.23b) 
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As  in  Equation  (3.16),  the  relationship  between  the  inward  and  outward  propagating 
waves  in  region  1  can  be  obtained  by  relating  ft\  and  g]x  to  c „  and  d'n  by  writing 
Equation  (3.23)  in  matrix  form  as 


ZTEiW\Kr»)  -  WIK' o)  )(^r„)  -  HW’iklrv) 

~Zn,-Z-HW\k'nra)  -  Hm(k'„ru)  Zn^H^'(k'Lru)  +  H^(k'LT„) 


s'„ 


•WM  -  Zte£J„ (*>o)  r„(k'L?0)  -  ZTE^MKr„) 

Jn(k'Rr„)  +  Zru^JMn)  -  ZT„±JUk'Lru) 


(3.24) 


Solving  Equation  (3.24)  for  /„  and  g'n,  as  was  done  for  Equation  (3.16),  results 
in  [V°]  matrix  elements  given  by 

=  ),r„)  +  !  (  Zte Z„, 

A,  V  l] 

+  -%-  \ZTEUk\,r„)HVKklr„)  -  Zl  ttJ^k)tTa)H^\k\ru)]  (3.25a) 

V\&n  L 
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(3.25b) 


yii  __  _  (i  _  zteztm\  2 j 

"  V  Vl  )  A„ 

Jn(*>  r0)//(2)'(*^o)  +  J;(^o)ff<2)(*Ar0)  ,  Zte^TM  \ 

An  {  +  1?  ) 

+  -^~  [ZrEMklr^iKrn)  -  ZT*,J'Jk'Lrn)H«' »'(fcfcr„)] 

where 


(3.25c) 


(3.25d) 


A„  =  [aW(*ir„)aW'(*i.r„)  +  Jri2>'(Kr0)J/m(klr0)]  (l  +  — 1 ~'j 

-  ^  [ZreffW(*J,ru)ffl,,(*l'*)  -  Zr»/afl'(^r„)a(»'(^r„)]  .  (3.26) 

A  perfect  magnetic  conducting  center  cylinder  is  the  special  case  of  the  surface 
impedance  boundary  condition  where  Zje  —  ztm  — *  oo. 
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3.2.3  Material  Center  Cylinder 

If  there  is  no  center  impenetrable  cylinder,  i.e.,  region  1  extends  to  the  origin,  then 
the  requirement  of  finite  fields  at  *he  origin  forces  /,{  =  g\  =  0,  i.e.,  no  outward 
propagating  waves  in  region  1.  To  use  the  same  formulation  as  in  Equation  (3.20), 
define  [F°]  by 

10  0  0 
0  10  0 
0  0  0  0 
0  0  0  0 

3.2.4  Exterior  Region 

By  defining  the  exterior  region  to  be  layer  M  - 1-1,  the  relationship  between  the 
coefficients  for  layers  1  and  M  +  1  is  given  by 

A„  AI2  A, 3  Xu 
A2|  A22  A23  a2, 

A31  A'32  A33  Asj 

A4 1  A42  A,  3  Xu 

where  the  matrix  [A"]  is  given  by 

[a]  =  [ym]  [yw-]]  •  •  •  [: y 2]  [r1]  [r°] .  (3.29) 


Note  that  Equation  (3.28)  applies  for  chiral  cylinders  with  or  without  a  surface 
impedance  center  cylinder  by  the  proper  choice  of  [K0].  Requiring  the  different 
representations  for  fields  in  the  exterior  region  to  be  equal,  i.e.,  E1  -f  E5  =  Efl,+1 
and  H1  -f  H*  =  Hfl/+1,  results  in  the  following  relationships  between  ( a„,b„ )  and 


fM  + 1 
J  n 


TM 

=  (c^1  -  <•+')  J.(M 

+  (/„"+'  -»"*')  ^!,(M  (3.30a) 
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7 1: 


J^k,p)+KHW'{k,p)  = 

(c"+'  +  <£,+l) 

+  (fit'*' +  d'+')  nm'M 

(3.30b) 

TE 

X(iw)+b,M2}M  = 

(c;,+1  +  <'+')  -MM 

+  (/„*’+l  +g"+,)W,M 

(3.30c) 

TM 

JW^)+anHW'{k,p)  = 

(c»,+,  _  jjft.)  j;(fcop) 

(3.30d) 

TM  TE 

where  the  and  indicate  terms  that  are  nonzero  for  TMz  and  TE z  incident 
fields,  respectively.  From  Equation  (3.30)  the  coefficients  c*,+l  and  d*l  +  '  are  given 

by 

1 


Jl+ 1  _ 


-d*,+l  =  ±- 
2 


(3.31a) 

(3.31b) 


where  the  plus  and  minus  signs  of  ±  are  used  for  TM^  and  TE^  incident  fields, 
respectively. 

Using  Equations  (3.28)  and  (3.31)  the  solution  for  the  coefficients  c]n  and  d'n  in 
layer  1  is 

i  A22  i  X\2 


C„  = 


2(XuX22  —  X2\  An) 

=  A21±A„ 

2(AnAr22  —  A2i  An) 


(3.32a) 

(3.32b) 


where  the  plus  and  minus  signs  of  ±  are  used  for  TMz  and  TE z  incident  fields, 
respectively.  From  Equation  (3.30)  the  coefficients  for  the  scattered  fields  are 


SI/  +  1  _  M  + 1 
9n 


an  =  /' 

*>n  =  f?+'+£  +  ' 


(3.33a) 

(3.33b) 
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These  coefficients  have  the  following  properties: 


(3.34a) 

(3.34b) 

(3.34c) 


a_n  =  ar 


b- „  =  b„ 

an(TEz )  =  bn(TMz). 


As  a  result  of  Equation  (3.34c),  the  cross-polarized  echo  widths  of  circular  cylinders 
for  TM^  and  TEZ  incident  plane  waves  are  identical. 

In  summary,  Equation  (3.33)  provides  the  unknown  coefficients  in  the  eigenfunc¬ 
tion  expansions  for  the  scattered  fields  given  by  Equation  (3.5),  which  can  be  written 
in  the  far  field  as 


Ez 
E * 


-J*oP 


£  ant** 
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7T 


E  b' 


ein4> 


(3.35a) 

(3.35b) 


using  the  large  argument  approximation  for  the  Hankel  function.  In  addition,  the 
internal  fields  are  found  using  Equation  (3.7),  where  the  coefficients  are  found  using 


'  c  ' 
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fn 

=  [: Ym~' ]  [rm-2]  •• 

M  M  M 
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0 

3.3  Numerical  Results 


(3.36) 


This  section  presents  numerical  results  for  TMz  and  TE^  scattering  by  several  dif¬ 
ferent  multilayer  circular  chiral  cylinders  in  free  space.  Comparison  is  made  with 
similar  multilayer  achiral  cylinders  in  each  case.  For  convenience,  all  data  is  at  a 
frequency  of  300  MHz.  The  chiral  parameters  used  in  this  section  obey  the  limit  set 
in  [102]  of  \(c\  <  yJTJJi. 
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3.3.1  TM^  Incident  Plane  Wave 

This  section  presents  the  bistatic  scattering  patterns  for  a  low-loss  two-layer  chiral 
cylinder  with  and  without  a  perfectly  conducting  center  cylinder  for  a  TMz  incident 
plane  wave  with  comparison  to  equivalent  achiral  cylinders.  Figure  3.2  shows  the 
co-polarized  and  cross-polarized  bistatic  echo  widths  for  a  TMz  plane  wave  [E':  = 
e-j* o*)  incident  upon  the  chiral  cylinder  with  a  perfect  electric  conducting  center 
cylinder.  The  outer  radius  is  1.25  free  space  wavelengths.  The  curve  marked  ‘Achiral 
(£r  =  0)’  is  the  bistatic  echo  width  for  an  identical  multilayer  cylinder  with  the 
chirality  admittance  set  to  zero  in  both  layers.  The  chirality  of  the  material  clearly 
caused  a  significant  change  in  the  co-polarized  scattered  fields  while  generating  a  low 
level  cross-polarized  scattered  field.  The  chirality  admittance  in  the  inner  and  outer 
layers  is  65%  and  33%,  respectively,  of  the  theoretical  maximum  of  |£r|  <  Je/fi  set 
in  [102]. 

Figure  3.3  shows  the  co-polarized  and  cross-polarized  bistatic  echo  widths  for 
chiral  and  achiral  cylinders  with  the  same  geometry  as  those  in  Figure  3.2  except 
that  the  conducting  center  cylinder  has  been  replaced  by  the  same  material  as  the 
inner  layer.  Note  that  the  co-polarized  and  cross-polarized  fields  are  of  comparable 
magnitudes,  indicating  that  the  chirality  of  the  cylinder  has  caused  a  significant 
rotation  in  the  polarization  of  the  scattered  fields.  In  addition,  the  bistatic  pattern  of 
the  achiral  cylinder  has  more  lobes  than  that  of  the  chiral  cylinder.  This  is  unusual 
because  the  chiral  cylinder  is  larger  electrically,  since  one  of  the  two  wavelengths 
in  chiral  media  is  always  smaller  than  that  of  the  corresponding  achiral  media; 
therefore,  its  pattern  should  have  more  lobes.  Figure  3.4  shows  this  point  more 
clearly  using  a  polar  plot  where,  because  of  symmetry,  the  upper  half  shows  the 
co-  and  cross-polarized  bistatic  echo  width  for  the  chiral  cylinder  and  the  lower  half 
shows  the  echo  width  for  the  achiral  cylinder. 
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Cnirai  Cylinder  with  PlC 
Frequency  =  300  MHz 


Figure  3.2:  The  co-polarized  and  cross-polarized  bistatic  echo  width  of  a  two-layer 
chiral  cylinder  with  a  perfect  electric  conducting  center  cylinder  for  a  TMz  incident 
plane  wave. 


Two-layer  Chiral  Cylinder 

Frequency  =  300  MHz 


Figure  3.3:  The  co-polarized  and  cross-polarized  bistatic  echo  width  of  a  two-layer 
chiral  cylinder  for  a  TM^  incident  plane  wave. 


3.3.2  TE 7  Incident  Plane  Wave 


Tins  section  presents  the  bistatic  scattering  patterns  for  a  high-loss  two-layer  chiral 
cylinder  with  and  without  a  perfectly  conducting  center  cylinder  for  a  TE^  incident 
plane  wave.  Internal  fields  and  approximate  surface  currents  on  the  center  cylinder 
are  presented  for  the  chiral  cylinder  with  the  perfectly  conducting  center  cylinder. 

Figures  3.5  and  3.6  are  similar  to  Figures  3.2  and  3.3,  except  that  they  are  for 
the  TE7  polarization  and  for  higher  loss  chiral  media.  In  these  figures  the  cross- 
polarized  fields  are  greater  than  the  co-polarized  fields  for  angles  between  20  and  70 
degrees,  indicating  that  the  chirality  of  the  cylinder  has  caused  a  significant  rotation 
in  the  polarizatioh  of  the  scattered  fields.  The  chirality  admittance  in  the  inner  and 
outer  layers  is  23%  and  80%,  respectively,  of  the  theoretical  maximum  of  |£,|  <  yjt/fi 
set  in  [102]. 

Figures  3.7  and  3.8  show  the  near  zone  electric  and  magnetic  fields  along  the 
positive  y  axis  for  the  chiral  and  achiral  cylinders  of  Figure  3.5  with  the  perfect 
electric  conducting  center  cylinder.  At  the  two  material  boundaries,  y  =  0.75  and 
y  =  1.25  meters,  Ey  and  Hy  are  discontinuous  as  required.  At  the  boundary  between 
free  space  and  the  outer  material  layer,  the  Ex  and  Hz  fields  of  the  achiral  cylinder 
exhibit  a  much  smaller  slope  discontinuity  than  the  fields  of  the  chiral  cylinder. 
Interestingly,  these  same  fields  at  the  inner  material  boundary,  y  =  0.75  meters, 
exhibit  almost  no  slope  discontinuity  in  both  cylinders.  However,  the  related  E2  and 
Ht  fields  exhibit  sharp  slope  discontinuities  at  both  boundaries  of  the  chiral  cylinder. 
In  Figure  3.7,  the  magnitude  of  HT  at  the  surface  of  the  conducting  cylinder  is  greater 
than  that  of  Hz.  This  means  that  at  the  point  (x  =  0,y  =  0.25m),  the  J2  surface 
current  on  the  conducting  cylinder  is  greater  than  the  Jx,  i.e.,  the  J $  surface  current. 
For  a  conducting  cylinder  coated  with  an  achiral  material  only  the  J 0  surface  current 
will  be  nonzero  for  a  TE^  incident  plane  wave  as  shown  in  Figure  3.8. 

Figure  3.9  shows  the  tangential  magnetic  fields  at  the  surface  of  the  conducting 
center  cylinder  for  the  chiral  and  achiral  cylinders  shown  in  Figure  3.5.  The  Jt 
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Echo  Width  (dB/m) 


Chiral  Cylinder  with  PLC 
Frequency  =  300  MHz 


Figure  3.5:  The  co-polarized  and  cross-polarized  bistatic  echo  width  of  a  two  layer 
chiral  cylinder  with  a  perfect  electric  conducting  center  cylinder  for  a  TE^  incident 
plane  wave. 
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Two— layer  Chiral  Cylinder 
Frequency  =  300  MHz 

yf 


Figure  3.6:  The  co-polarized  and  cross-polarized  bistatic  echo  width  of  a  two  layer 
chiral  cylinder  for  a  TEz  incident  plane  wave. 


7^|H|  (V/m)  |E|  (v/m) 


Figure  3.7:  Internal  fields  along  x  =  0  for  the  chiral  cylinder  shown  in  Figure  3.5. 
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current  is  greater  than  the  J $  current  over  the  entire  surface  of  the  conducting 
cylinder  in  the  center  of  the  two-layer  chiral  cylinder,  even  though  the  J2  current 
is  zero  on  the  surface  of  a  conducting  cylinder  at  the  center  of  an  identical  achiral 
cylinder. 

Figure  3.10  shows  the  near  zone  electric  and  magnetic  fields  along  the  x  axis, 
for  the  geometry  shown  at  the  top  of  Figure  3.5  with  the  perfect  electric  conducting 
center  cylinder. 

3.3.3  Chirality  Versus  Scattering 

This  section  examines  the  effect  of  the  chirality  admittance  on  the  backscatter  and 
bistatic  scattering  patterns.  First,  a  homogeneous  chiral  cylinder  with  fi,  =  2,  er  =  3, 
and  a  radius  of  0.15  meters  is  examined  over  the  range  of  £c  =  0  (achiral)  to  the 
theoretical  maximum  of  161  <  =  0.00325  [102],  for  lossless  and  lossy  media. 

Figure  3.11  and  3.12  present  the  co-  and  cross-polarized  backscatter  echo  width 
versus  the  magnitude  of  the  chirality  admittance  for  TMz  and  TEz  incident  plane 
waves.  In  all  four  cases,  the  cross-polarized  echo  width  increases  rapidly  as  (r 
increases  and  equals  the  co-polarized  echo  width  when  £c  is  approximately  18%  of 
the  theoretical  maximum.  The  most  interesting  case  of  the  four  is  the  TEz  plane 
wave  incident  on  the  loss/  chiral  cylinder.  In  this  problem  the  backscattered  field 
is  entirely  cross-polarized  for  a  small  range  of  the  chirality  admittance.  Figure  3.13 
show's  the  co-polarized  bistatic  scattering  pattern  of  this  cylinder  for  a  TEz  incident 
plane  wave  for  the  full  range  of  £c.  Even  through  the  achiral  (£r  =  0)  bistatic 
scattering  pattern  has  no  deep  nulls,  this  figure  shows  that  deep  nulls  will  occur 
for  certain  ranges  of  6-  Interestingly,  nothing  significant  happens  at  the  theoretical 
maximum  of  |£r]  <  yjTffi  set  in  [102]. 

Next,  a  homogeneous  chiral  cylinder  with  fir  =  2,  er  =  3,  a  radius  of  0.15  meters, 
and  a  perfectly  conducting  center  cylinder  of  radius  0.05  meters  is  examined  over  the 
range  of  £c.  Figures  3.14  and  3.15  present  the  co-  and  cross- polarized  backscatter 
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Chiral  Cylinder  with  PEC 


Frequency  =  300  MHz 


Figure  3.10:  Internal  fields  along  y  =  0  for  the  chiral  cylinder  shown  in  Figure  3.5 
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Backscatter  versus  Chirality 
er  =  3,  jUr  =  2,  radius  =  0.1  5  m 
Frequency  =  300  MHz 


Figure  3.11:  Backscatter  echo  width  versus  the  magnitude  of  the  chirality  admittance 
for  TMz  incident  plane  wave. 


Echo 


Figure  3.12:  Backscatter  echo  width  versus  the  magnitude  of  the  chirality  admittance 
for  TE z  incident  plane  wave. 


B i stall c  Scattering  versus  Chirality 


IJLr=  2,  tan  <5m=0.05 
er=3,  tan  <5e=0.05 
Radius  =  0.15  m 


Freq.  =  300  MHz 


0.001 


Figure  3.13:  Bistatic  echo  width  versus  the  magnitude  of  the  chirality  admittance 
for  TEz  incident  plane  wave. 
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echo  width  versus  the  magnitude  of  the  chirality  admittance  for  TMz  and  TE z  inci¬ 
dent  plane  waves.  In  all  four  cases,  the  cross- polarized  echo  width  increases  rapidly 
as  £r  increases  and  equals  the  co-polarized  echo  width  when  £c  is  approximately  40% 
of  the  theoretical  maximum.  As  the  radius  of  the  conducting  cylinder  is  increased, 
the  thickness  of  the  coating  is  reduced  and  less  material  interacts  with  the  incident 
field,  thereby  producing  less  cross-polarized  scattered  fields.  This  relationship  is 
easily  seen  when  comparing  Figures  3.11  and  3.12  with  Figures  3.14  and  3.15. 
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3.4  Summary 


This  chapter  has  presented  an  efficient  recursive  eigenfunction  solution  for  scattering 
by  a  general  multilayer  circular  chiral  cylinder,  with  or  without  a  surface  impedance 
center  cylinder,  for  a  TM^  or  TE z  normal  incidence  plane  wave.  The  internal  fields 
were  expanded  as  inward  and  outward  propagating  right  and  left  circularly  polarized 
waves,  using  circularly  polarized  vector  wave  functions,  resulting  in  four  unknown 
coefficients  in  each  layer.  The  resulting  solution  requires  the  multiplication  of  M  4x4 
matrices  for  an  M  layer  cylinder.  The  TE z  solution  differs  from  the  TM?  solution 
by  only  a  plus/minus  sign  at  one  step.  Including  a  surface  impedance  center  cylinder 
requires  only  one  additional  matrix  multiplication.  An  examination  of  coefficients 
responsible  for  the  scattered  fields  showed  that  the  cross-polarized  bistatic  echo 
width  of  circular  cylinders  is  identical  for  TM^  and  TE;?  incident  plane  waves. 

TM;?  and  TE;?  bistatic  scattering  patterns  for  two-layer  chiral  cylinders  with  and 
without  a  perfectly  conducting  center  cylinder  were  presented  as  reference  solutions. 
These  patterns  demonstrated  that  the  chirality  of  these  cylinders  caused  a  significant 
rotation  in  the  polarization  of  the  scattered  fields.  More  significant  was  the  effect 
of  the  chirality  on  the  internal  fields.  In  one  case,  more  J2  than  J &  current  was 
induced  on  the  perfectly  conducting  cylinder  coated  by  two  layers  of  chiral  media 
for  a  TEz  incident  plane  wave,  even  though  the  achiral  problem  would  have  induced 
no  Jz  current  on  the  conducting  center  cylinder.  The  chirality  admittance  in  these 
cases  was  greater  than  50%  of  the  theoretical  maximum  of  [£c(  <  \Je/n  set  in  [102]. 

To  examine  the  relationship  between  the  magnitude  of  the  chirality  admittance 
and  scattered  fields,  backscatter  and  bistatic  scattering  patterns  were  presented  for 
the  full  range  of  ( c .  Even  through  the  achiral  bistatic  scattering  pattern  may  have 
no  deep  nulls,  the  co-polarized  bistatic  pattern  for  a  chiral  cylinder  can  have  deep 
nulls  for  certain  ranges  of  (c.  In  one  case  the  backscattered  field  was  entirely  cross- 
polarized  for  a  small  range  of  the  chirality  admittance.  The  cross-polarized  backscat¬ 
tered  fields  increased  in  similar  fashion  for  all  of  the  geometries  examined  as  £r  was 
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increased.  1  ne  rale  of  increase  was  dependent  on  the  amount  of  chiral  material, 

since  chirality  is  a  bulk  effect. 
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Chapter  4 

Scattering  by  an  Arbitrary 
Cylinder 


This  chapter  presents  a  pulse-basis  point-matching  method  of  moments  solution  [123, 
124]  to  the  two-dimensional  problem  of  scattering  by  an  isotropic  inhomogeneous 
chiral  cylinder  of  arbitrary  cross  section  in  free  space  or  in  the  presence  of  a  perfectly 
conducting  half-plane.  The  method  of  moments  solution  for  scattering  by  a  chiral 
cylinder  is  obtained  using  the  volume  equivalence  theorem  developed  in  Section  2.4  to 
replace  the  chiral  scatterer  by  free  space  and  equivalent  electric  and  magnetic  volume 
polarization  currents  (J,  M),  as  shown  in  Figure  4.1.  By  enforcing  the  volume 
equivalence  theorem  in  the  chiral  scatterer,  a  pair  of  coupled  vector  integral  equations 
for  the  currents  (J,  M)  are  obtained.  These  r-,upled  vector  integral  equations  are 
equivalent  to  six  coupled  scalar  integral  equ;  iio-  for  ( Jz ,  Jy,  Jz)  and  (A/r,  My,  A/-), 
which  are  solved  using  a  pulse-basis  point-matching  method  of  moments  solution. 
This  problem  could  be  formulated  in  terms  of  surface  currents  on  the  chiral  cylinder 
[74,  80,  95].  The  main  advantage  of  the  volume  formulation  is  that  it  can  more  easily 
Lr  »  inhomogeneous  media. 

The  electric  surface  currents  on  the  half-plane  are  not  explicitly  included  as 
unknowns  in  the  method  of  moments  solution.  Instead,  their  effects  are  exactly 
accounted  for  by  including  the  half-plane  Green’s  function  in  the  kernel  of  the  integral 
equation  [125,  sec.  8.3].  Efficient  techniques  for  accurately  evaluating  the  integrals 
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Figure  4.1:  Volume  Equivalence  Theorem  (a)  Original  problem:  impressed  currents, 
half-plane,  and  chiral  scatterer.  (b)  Equivalent  problem:  scatterer  replaced  by  free 
space  and  the  equivalent  currents. 

in  this  Green’s  function  are  presented  in  Appendix  B.  In  the  free  space  problem  the 
half-plane  Green’s  f  jnction  is  replaced  by  the  free  space  Green’s  function. 

4.1  Derivation  of  General  Integral  Equations 

This  section  develops  a  set  of  coupled  integral  equations  for  the  equivalent  electric 
and  magnetic  currents  representing  a  chiral  scatterer  in  free  space  or  in  the  presence 
of  a  perfectly  conducting  half-plane.  In  the  original  problem  of  Figure  4.1(a),  the 
impressed  currents  (J1,  M1)  radiate  the  total  fields  (E,  H)  in  a  medium  which  is 
free  space  except  for  the  perfectly  conducting  half-plane  and  a  chiral  scatterer  with 
constitutive  parameters  (/x,  e,  £c)  confined  to  the  region  R.  The  scatterer  may 
be  lossy  and  inhomogeneous.  As  shown  in  Figure  4.1(b),  the  volume  equivalence 
theorem  developed  in  Section  2.4  is  used  to  replace  the  chiral  scatterer  by  free  space 
and  the  equivalent  electric  and  magnetic  volume  polarization  currents 

J  =  yw(ec  -  e0)E  4-  H  (4.1a) 

M  =  —  /xn)H  —  u;/x£rE  (4-lb) 


confined  to  the  region  R. 


57 


In  the  equivalent  problem  of  Figure  4.1(b),  the  total  fields  (E,  H)  at  any  point 
in  space  are  given  by 

E  =  E'  +  EJ  +  Em  (4.2a) 

H  =  H'  +  Hj  +  Hm  (4.2b) 


where  (E\  H‘),  (EJ,  HJ),  and  (EM,  HM)  are  the  fields  radiated  by  (J1,  M'),  J,  and 
M,  respectively,  in  free  space  or  in  the  presence  of  the  half- plane.  Then,  substituting 
Equation  (4.2)  into  Equation  (4.1)  yields 


-  EJ  -  EM  -  c£  [HJ  +  Hm]  +  aEJ  = 

E'  +  ceH' 

(4.3a) 

-H3  -  HM  -  cm  [E3  +  EM]  +  aMM  = 

H’  +  cmE’ 

(4.3b) 

in  region  R,  where  cE,  cm,  aE,  and  am  are  defined  by 

“m(c 

ju>(ec  -  eu) 

(4.4a) 

MM  ~  Mo) 

(4.4b) 

1 

M€c  ~  e„) 

(4.4c) 

1 

aM  Mm -Mo)' 

(4.4d) 

Equation  (4.3)  can  be  considered  to  be  a  coupled  integral  equation  for  (J,  M)  since 
the  fields  (EJ,  HJ)  or  (EM,  HM)  can  be  written  as  integrals,  over  the  region  R 
containing  the  chiral  scatter,  of  J  or  M,  respectively,  dotted  into  the  appropriate 
dyadic  Green’s  function  [126].  For  example,  EM  could  be  written  as 

EM(r)  =  j  M(r')  •  G(r,  r')  dv ,  (4.5) 

R 

where  r'  is  the  source  point,  r  is  the  field  point,  and  G  is  the  electric  field  dyadic 
Green’s  function  for  the  magnetic  source  radiating  in  free  space  or  in  the  presence 
of  a  perfectly  conducting  half-plane. 
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4.2  Integral  Equations  for  a  Two-Dimensional 
Body 

In  this  section,  the  general  integral  equation  of  the  previous  section,  Equation  (4.3), 
is  simplified  for  the  case  of  normal  incidence  scattering  by  chiral  cylinders  of  arbitrary 
cross  section.  Since  this  is  now  a  two-dimensional  problem,  there  are  some  simpli¬ 
fications  from  the  general  three-dimensional  case.  In  particular,  the  TM  currents 
(Jz,  Mx,  My)  generate  only  the  ( Ez ,  Hx,  Hy )  fields  and  the  TE  currents  (A/*,  Jx,  Jv) 
generate  only  the  ( Hz ,  Ex,  Ey)  fields.  In  this  case,  Equation  (4.3)  reduces  to  the 
following  six  scalar  equations  valid  in  R: 

-  E'z  -  E"*  -  E'z'v  -  ce(H*u  +  HJZ '  +  HJZ«)  +  aEJz  =  E\  4-  cEH\  (4.6a) 

- Hi '  -  Hilz  -  HAT’«  -  cm{E*'‘  +  El *  +  EJr »)  4-  aMMx  =  H'x  +  cME'x  (4.6b) 

-Hl‘  -  H"r  -  H*lv  -  cm{E*'z  +  Ef  4-  El')  4-  aMMv  =  H\,  +  chlE[  (4.6c) 

-cm{Ej/  +  E?1'  +  Ei'v)  -  Hi'*  -  HJZ*  -  Hi*  +  aMM2  =  H\  4-  cmE[  (4.6d) 

-cE{Hl*  4-  Hi'*  +  Hi'*)  -  Ei,z  -  Ei*  -  EJX*  4-  aEJx  =  E'x  +  cEH'x  (4.6e) 

-cE(Hl*  +  Hi1*  +  Hi’>)  -  Ei,z  -  EJy*  -  EJy*  +  aEJy  =  E'y  4-  cEH'y.  (4.6f) 

In  Equation  (4.6),  the  superscript  indicates  the  source  and  the  subscript  the  field 
component.  For  example,  H l*  is  the  y  component  of  the  magnetic  field  of  Jz  and  E'x 
is  the  x  component  of  the  electric  field  of  (J\  M1)  in  free  space  or  in  the  presence 
of  a  half-plane.  For  normal  incidence  scattering  by  achiral  cylinders,  the  TM  and 
TE  currents  are  uncoupled  and  can  be  determined  separately  [127]— [129] .  However, 
Equation  (4.6)  shows  that  the  TM  and  TE  currents  are  coupled  for  chiral  cylinders, 
and  it  is  therefore  necessary  to  solve  for  all  currents  simultaneously. 
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Figure  4.2:  A  trapezoidal  cross  section  cylinder  split  into  N  smaller  trapezoidal  cells. 


4.3  Method  of  Moments  Solution 


The  pulse-basis  point-matching  method  of  moments  solution  of  Equation  (4.6)  is 
presented  in  this  section.  For  convenience,  the  chiral  cylinder  is  restricted  to  be 
composed  of  one  or  more  blocks  of  trapezoidal  cross  section.  Without  loss  of  gener¬ 
ality,  consider  a  single  trapezoidal  cross  section  cylinder  as  shown  in  Figure  4.2.  The 
trapezoid  is  split  into  N  smaller  trapezoidal  cells  by  dividing  the  opposite  sides  of 
the  trapezoid  into  equal  segments.  These  cells  must  be  sufficiently  small  so  that  the 
electromagnetic  fields  and  the  constitutive  parameters  of  the  chiral  cylinder  can  be 
considered  as  essentially  constant  within  each  cell.  For  a  pulse-basis  point-matching 
method  of  moments  solution  of  scattering  by  achiral  cylinders,  the  maximum  dimen¬ 
sion  of  the  cells  is  typically  less  than  one-tenth  of  the  wavelength  in  the  material. 
The  region  occupied  by  the  nth  cell  is  denoted  by  R„ ,  and  has  cross-sectional  area 
An.  Now  define  the  volume  electric  and  magnetic  current  pulse  expansion  functions 
for  n  =  1,2, . . . ,  N  as 


Pi 


1  /j4Ti  A/m2,  within  R„ 

< 

0,  otherwise 


(4.7a) 


< 


within  Rn 
otherwise. 


(4.7b) 


pM 

r  rx 


1  /A„  V/nT, 

0, 


The  equivalent  electric  and  magnetic  currents  can  then  be  expanded  as 


J  =  xjx  +  yJv  +  zJ2 

=  i'£c^  +  y'tcy-pn+i'tc„p^  (4.8a) 

n=  1  n=  1  n=  I 

M  =  xMx  -f  yMv  -f  zM2 

=  x-£K,„P"  +  y'£KmPZ'  +  z-£K,„P!'.  (4.8b) 

n= I  n— 1  n=  J 

Note  that  for  N  cells  there  are  QN  unknown  coefficients:  Cxn ,  Cyn ,  Czn ,  Kxn ,  Kyn, 
and  K2n  (n  =  1,2,...,  N). 

The  pulse-basis  point-matching  method  of  moments  solution  is  obtained  by  sub¬ 
stituting  Equation  (4.8)  into  Equation  (4.6)  and  requiring  the  resulting  six  equations 
to  be  exactly  satisfied  at  the  centroid  of  each  of  the  N  cells.  This  transforms  Equa¬ 
tion  (4.6)  into  a  6 TV  x  6 N  system  of  simultaneous  linear  equations  which  can  be 
compactly  written  in  matrix  form  as 


[Z  +  AZ]I  =  V. 


(4.9) 


This  equation  cm  be  symbolically  decomposed  into  submatrices  representing  the 
relationships  between  the  TMz  and  TE^  fields  and  currents 


TM/TM 

TM/TE 

TM 

TM 

TE/TM 

TE/TE 

TE 

TE 

(4.10) 
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4.3.1  The  Impedance  Matrix  [Z\ 

For  an  inhomogeneous  chiral  cylinder,  the  impedance  matrix  \Z\  can  be  written  as 
a  6  x  6  block  matrix 


(Z]  = 


Zu 

Z  i  2 

Z\3 

CeZ41 

CeZ15 

CeZxg 

Z21 

Z22 

Z23 

CmZ^x 

CmZss 

Cm  Z$g 

z,  1 

Z32 

Z33 

Cm  Z&4 

CmZg; 

CmZg6 

CmZ\\ 

CmZ\2 

CmZi3 

Z.U 

Zx  3 

Z  XG 

Cr.Z2\ 

CEZ22 

CEZ23 

Z54 

Z55 

Zgg 

CEZM 

CeZ'32 

CEZs3 

Zq4 

Z$5 

Zg6 

(4.11) 


where  each  of  the  ZtJ  is  an  N  X  N  matrix,  where  the  m,n  element  is  the  negative 
F,  field  of  the  P_,(n)  current  expansion  function  evaluated  at  the  centroid  of  the  cell 
m,  i.e., 

~  ~ •f|'(P,;(n))|m  (4-12) 

where 


F,  =  { 


Ez, 

i 

=  1 

Hr, 

i 

=  2 

Hy, 

i 

=  3 

Hz, 

i 

=  4 

Ex, 

i 

-  5 

Ey , 

i 

=  6 

(4.13) 
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and 


i(") 


=  < 


(4.14) 


*Pn,  J  =  1 
iP",  j  =  2 
yPn,  j  =  3 
J=4 

*Pt,  J  =  5 
l  yPS,  j  =  6 

for  m,n  =  1, 2, . . . ,  iV.  For  example,  the  m,n  element  of  Z23  is  —Hx(yP^')  evaluated 
at  the  centroid  of  cell  m. 

The  Ce  and  Cm  are  N  x  N  diagonal  matrices  given  by 

0 


Ce  = 


Ce\ 


CE 


A’  J 


ii; 
(i.iO 


a) 


Ca/  = 


cm 


cm 


k  J 


(4.15b) 


where  ce\  is  ce  from  Equation  (4.4)  evaluated  at  the  centroid  of  cell  m  and  similarly 

I  m 

for  cm  . 

m 

The  computational  effort  to  obtain  the  free  space  impedance  matrix  [Z]  can  be 
reduced  by  recognizing  simple  relationships  between  the  block  matrices.  Let 
(EJ,HJ)  and  (EM,HM)  denote  the  free  space  fields  of  the  arbitrary  electric  and 
magnetic, currents  J  and  M,  respectively.  Then,  directly  from  Maxwell’s  equations, 
it  can  be  shown  that  if  J  is  numerically  equal  to  M,  then 

EM  =  -HJ  (4.16a) 

HM  =  EJ/V20,  (4.16b) 

where  ?7U  =  is  the  characteristic  impedance  of  free  space  and  =  should  be 

read  “numerically  equal”.  Examining  Equations  (4.13)  and  (4.14),  and  noting  that 
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Fr!  =  P",  yields 


"  1 "  1 

Z45 

Zag 

iZu 

—  Z12 

1  . 

1 

Z54 

Z55 

Z56 

= 

— Z21 

Vi,Z22 

Vu  Z24 

^61 

^65 

Z&6 

nlz,. 

4.3.2  The  Impedance  Matrix  [A Z] 


The  impedance  matrix  [A Z]  is  a  67V  x  67V  diagonal  matrix  given  in  block  matrix 
form  as 

r  AZn  0 

A  Z22  0 

0  A  Z33 


[AZ]  = 


A  Zu  0 

0  I  AZ55 

0  A  Zee 

Each  of  the  AZ„’s  is  a  TV  x  TV  diagonal  matrix  given  by 


AZ„  = 


a,  JA\ 


0 


ai  .,M A' 

A 


where  Am  is  the  cross-sectional  area  of  cell  m  and  a,  is  given  by 


a, 


0£-  ,  for  t=l,5,6 

Irn 

Ofl/1  ,  for  i=2, 3, 4, 


(4.18) 


(4.19) 


(4.20) 


where  a  A  is  a  1;  from  Equation  (4.4)  evaluated  at  the  centroid  of  cell  m  and  similarly 

'\rn  ' 

for  d/i/l  . 
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4.3.3  The  Current  Vector  I 


The  current  vector  I  is  defined  as  a  6  element  block  vector 


cz  ' 

Kx 


Kz 
Cx 
.  cv 


where  the  coefficients  Cxn ,  Cyn,  Czn,  Kxn,  Kyn ,  and  Kzn  (n  = 
tion  (4.8)  are  contained  in  the  length  N  column  vectors  Cx , 
Kz ,  respectively. 


(4.21) 


1,2, ... ,  N )  of  Equa- 
Cy,  C„  Kx,  Kv,  and 


4.3.4  The  Voltage  Vector  V 


The  voltage  vector  V  is  also  defined  as  a  6  element  block  vector 


E\  +  CeH\ 
H'x  +  CmE'x 
H'y  +  CME'y 
H'z  +  Cm  E'z 
E'x  +  CeH'x 

E'y  +  CEH'y 


(4.22) 


where  the  terms  E'x,  E'y,  E'z,  H'x,  H'v,  H'z  are  length  N  column  vectors  consisting  of 
the  fields  of  the  impressed  currents  (J',M')  evaluated  at  the  centroid  of  cell  m  = 
1,...,7V.  For  example,  the  mth  element  of  E'z  is  the  z  component  of  the  electric 
field  generated  by  the  currents  (J',M'),  evaluated  at  the  centroid  of  cell  m.  CE  and 
Cm  are  the  N  x  N  diagonal  matrices  defined  in  Equation  (4.15). 
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y 


?e' 


}— - 1.2  m - 0.82  m  - 1.2  m 

jur  =  4.0,  er  =  1.5,  £c  =  -0.0016 

maximum  thickness  =  0.04  m 

Figure  4.3:  A  TM/  plane  wave  normally  incident  upon  a  tapered  chiral  slab. 

4.3.5  Scattered  and  Internal  Fields 

Equation  (4.9)  can  be  solved  for  the  current  vector  /  which,  when  substituted  into 
Equation  (4.8),  provides  an  approximation  to  the  equivalent  electric  and  magnetic 
currents  (J,M).  The  scattered  fields  are  simply  the  fields  of  (J,M). 

The  total  internal  fields  at  the  centroid  of  a  cell  are  given  by 

E  =  -J(M  ~  /*o)J  +  M&M 

«  [(«c  -  c0)(m  -  Po)  ~  P2tc2] 

„  -  €„)M  +  fi(c3 

H  -  r  T 

W  [(ec  -  Co)(M  -  Mu)  ~  P2tc2\ 
where  (p,ec,£c)  and  (J,M)  are  evaluated  at  the  centroid  of  the  cell. 

4.4  Free-Space  Numerical  Results 

This  section  presents  numerical  results  for  TM  and  TE  scattering  by  a  variety  of 
chiral  cylinders  in  free  space.  All  data  is  at  a  frequency  of  300  MHz,  and  includes 
both  echo  width  and  internal  fields. 

Figure  4.3  shows  a  TM/  polarized  plane  wave  normally  incident  upon  a  chiral  slab 
(£'  =  ejkoy,  where  fc„  is  the  free  space  wave  number).  The  center  section  of  the  slab 


(4.23a) 

(4.23b) 
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is  0.82  meters  wide  and  (j.Ul  meters  thui;.  The  lapered  sections  are  each  1.2  meters 

long.  The  slab  is  lossless  with  parameters  of  er  =  1.5,  fir  =  4.0,  and  £c  =  —0.0016 
A/V.  Near  the  origin  the  fields  for  the  tapered  slab  should  approximate  that  of  an 
infinite  slab  of  uniform  thickness  0.04  m.  Figures  4.4  and  4.5  show  a  comparison  of 
the  method  of  moments  solution  for  the  total  fields  of  the  tapered  slab  to  the  exact 
solution  developed  in  Appendix  C  for  the  infinite  uniform  slab  along  the  line  x  =  0. 
The  phase  of  the  fields  has  been  normalized  by  dividing  fields  by  e^oy,  i.e.,  the  phase 
of  a  plane  wave  propagating  in  free  space  in  the  —  y  direction.  The  incident  electric 
field  is  a  i  polarized  plane  wave.  When  this  wave  hits  the  uniform  slab  it  produces 
a  purely  z  polarized  reflected  wave,  resulting  in  a  partial  standing  wave  pattern  to 
the  right  of  the  slab.  As  the  wave  enters  the  slab,  the  polarization  begins  to  rotate, 
causing  the  sc  component  of  the  electric  field  to  increase,  and  the  z  component  to 
decrease.  Finally,  as  the  wave  exits  the  slab  it  is  a  linearly  polarized  wave,  but  with 
its  polarization  rotated  from  z.  The  solution  for  the  tapered  slab  has  a  small  sc 
polarized  back  scattered  field  which  is  not  present  for  an  infinite  slab. 

Next  the  method  of  moments  solution  is  compared  to  the  exact  eigenfunction 
solution  [20]  for  the  bistatic  echo  width  of  a  0.1  m  radius  circular  chiral  cylinder. 
The  inset  in  Figure  4.6  shows  a  TM^  plane  wave  (E\  =  e~’k°T)  incident  upon  a  cir¬ 
cular  cylinder  approximated  by  a  number  of  trapezoidal  cells,  corresponding  to  the 
method  of  moments  expansion  functions.  Denoting  d  as  the  maximum  dimension  of 
the  cells,  Figure  4.6  shows  the  eigenfunction  and  method  of  moments  solution  for 
d  ~  0.2Am  (18  cells)  and  d  —  0.lAm  (63  cells).  Am  denotes  the  minimum  wavelength 
for  left  or  right  circularly  polarized  plane  waves  in  the  chiral  media.  For  the  ma¬ 
terial  parameters  in  Figure  4.6,  Am  =  0.301  meters.  Note  that  as  the  number  of 
cells  increases,  the  method  of  moments  solution  approaches  the  exact  eigenfunction 
solution  for  both  the  co-polarized  and  cross-polarized  scattered  fields. 

A  comparison  of  the  eigenfunction  and  method  of  moments  solution  for  the 
bistatic  echo  width  for  a  TE?  plane  wave  ( H\  =  e~^T  j 77,,)  incident  upon  a  0.15 
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Tapered  Chiral  Slab 
jur  =  4.0,  er  =  1.5,  (c  = -0.00 16 
Freq.  =  300  MHz 


- Infinite  Slab 

- Moment  Method 


Chiral  Cylinder,  Radius  =  0.1  meters 
/xr  =  4.0,  er=  1.5,  £o=0.0005 
Freq.  =  300  MHz 


4>  (deg) 

Figure  4.6:  The  co-polarized  and  cross-polarized  bistatic  echo  width  of  a  circular 
chiral  cylinder  for  a  TM^  incident  plane  wave. 
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Chiral  Cylinder,  radius  =  0.1  5  meters 
/ir  =  2,  tan  <5m=  .05,  Freq.  =  300  MHz 
6r  =  3,  tan  de=.05,  £c=.002 


Figure  4.7:  The  co-polarized  and  cross-polarized  bistatic  echo  width  of  a  circular 
chiral  cylinder  for  a  TE?  incident  plane  wave. 


m  radius  lossy  circular  chiral  cylinder  is  shown  in  Figure  4.7.  Note  that  the  cross- 
polarized  field  is  greater  for  most  angles  than  the  co-polarized  field,  indicating  that 
the  chirality  of  the  cylinder  has  caused  a  significant  rotation  in  the  polarization  of 
the  scattered  fields.  For  comparison,  the  echo  width  for  an  achiral  cylinder,  obtained 
by  setting  £c  —  0,  is  also  shown.  A  segment  size  of  0.05  Am  was  required  to  obtain  a 
well  converged  result,  where  Am  =  0.228  meters  for  the  material  parameters  shown. 
The  backscatter  ( <f>  =  180°)  echo  width  for  this  geometry  was  highly  sensitive  to  the 
segment  size  and  in  Figure  4.7  the  greatest  error  occurs  at  <j>  =  180°.  The  backscat¬ 
ter  echo  width  of  this  cylinder  was  examined  in  Figure  3.12  of  Chapter  3  for  the 
full  range  of  and  found  to  be  highly  sensitive  to  the  value  of  (c  in  the  region  of 
£r  =  0.002.  The  sensitivity  of  this  cylinder  to  segment  size  may  be  related  to  the 
sensitivity  in  in  the  same  way  the  calculated  scattering  from  a  resonant  body  is 
sensitive  to  minor  changes  in  its  materials  and  geometry,  as  well  as  the  method  used 
to  calculate  the  scattering. 

The  magnitude  of  the  internal  fields  along  the  cylinder’s  centerline  are  shown  in 
Figure  4.8.  As  the  incident  Hz  field  propagates  through  the  cylinder,  the  Hz  internal 
field  decreases  and  the  Ez  field  increases  due  to  the  chirality  of  the  material.  This  is 
in  addition  to  the  effects  due  to  the  cylindrical  geometry.  The  most  significant  effect 
of  the  chirality,  besides  generating  cross-polarized  fields,  is  to  reduce  the  effective 
wavelength  in  the  cylinder,  as  well  as  the  actual  wavelength  in  the  medium.  For 
the  achiral  cylinder,  the  actual  material  wavelength  is  0.408  meters  with  a  spacing 
of  0.12  to  0.15  between  the  peaks  and  valleys  of  the  internal  fields.  The  right  and 
left  wavelengths  of  the  chiral  medium  are  0.228  and  0.730  meters,  respectively.  The 
spacing  between  peaks  and  valleys  of  the  internal  fields  varied  between  0.03  and  0.09 
meters,  with  0.05  meters  being  typical.  Based  on  the  achiral  cylinder,  the  spacing 
should  have  been  0.067  to  0.084  meters.  The  chirality  caused  approximately  a  25% 
reduction  in  the  effective  wavelength  in  the  medium  beyond  that  predicted  based 
on  the  smallest  chiral  wavelength.  This  becomes  very  important  when  choosing  the 
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Figure  4.8:  Internal  fields  along  y  =  0  for  the  chiral  cylinder  shown  in  Figure  4.7. 


cell  size  to  use  for  a  method  of  moments  solution.  For  the  point-matching  pulse- 
basis  solution  used  here,  the  cell  size  for  the  achiral  cylinder  should  have  been  0.041 
meters,  using  a  ten  cells  per  wavelength  rule.  Using  the  same  rule,  the  cell  size  in 
the  chiral  cylinder  should  have  been  0.023  meters.  However,  after  adjusting  for  the 
observed  reduction  in  the  effective  wavelength  in  the  chiral  cylinder,  the  cell  size 
should  be  0.017  meters.  Overall,  the  chirality  caused  a  factor  of  5  increase  in  the 
number  of  cells  needed  for  the  method  of  moments  solution. 

Figure  4.9  shows  the  bistatic  echo  width  for  a  0.1  m  x  2.0  m  chiral  slab  with  a 
TMz  plane  wave  incident  at  30°  from  the  x  —  axis.  The  cell  size  of  the  basis  functions 
is  0.057  Am,  where  A„,  =  0.35  meters  for  the  material  parameters  shown.  As  in  the 
previous  figures,  the  chirality  has  generated  a  substantial  cross-polarized  field.  The 
magnitude  of  the  internal  fields  along  the  sc  — axis  of  the  slab  are  shown  in  Figure  4.10. 
As  the  E2  incident  field  propagates  along  the  slab,  a  Hz  internal  field  is  generated 
by  the  chirality  of  the  material.  Figures  4.11  and  4.12  show  the  bistatic  echo  width 
and  internal  fields  for  the  0.1  m  x  2.0  m  chiral  slab  shown  in  Figure  4.9  for  a  TE^ 
plane  wave  incident  from  30°  off  the  x— axis.  The  co-polarized  bistatic  pattern  in 
Figure  4.11  is  very  similar  to  the  co-polarized  bistatic  pattern  in  Figure  4.9.  The 
cross-polarized  bistatic  patterns  are  almost  identical.  In  Figure  4.12,  the  internal 
fields  also  show  almost  the  same  patterns  as  in  Figure  4.10,  except  E  and  H  are 
reversed.  As  the  Hz  incident  field  propagates  along  the  slab,  a  Ez  internal  field  is 
generated  by  the  chirality  of  the  material,  just  as  in  Figure  4.10. 

Figure  4.13  shows  an  inhomogeneous  chiral  slab,  identical  to  that  shown  in  Fig¬ 
ure  4.9  except  that  (r  =  0.0005  A/V  for  x  >  0,  and  (c  =  —0.0005  A/V  for  x  <  0, 
with  a  TM^  incident  plane  wave  from  30°  off  the  x-axis.  The  most  significant 
change  from  Figure  4.9  is  that  the  field  scattered  back  along  the  x  — axis  is  20  dB 
lower.  The  magnitude  of  the  internal  fields  along  the  x-axis  of  the  slab  is  shown  in 
Figure  4.14,  which  is  similar  to  Figure  4.10  except  that  the  Hz  field,  which  increased 
along  the  slab  in  Figure  4.10,  decreases  suddenly  at  the  interface  between  the  two 
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Figure  4.9:  The  co-polarized  and  cross-polarized  bistatic  echo  width  of  a  0.1  m  x  2.0 
m  homogeneous  chiral  slab  for  a  TM?  plane  wave  incident  at  30°  off  the  x-axis. 
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Figure  4.11:  The  co-polarized  and  cross-polarized  bistatic  echo  width  of  a  0.1  m  x 
2.0  m  homogeneous  chiral  slab  for  a  TE^  plane  wave  incident  at  30°  off  the  x-axis. 
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Figure  4.12:  Internal  fields  along  y  =  0  for  the  homogeneous  chiral  slab  shown  in 
Figure  4.11. 
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Figure  4.13:  The  co-polarized  and  cross-polarized  bistatic  echo  width  of  a  0.1  m  x 
2.0  m  inhomogeneous  chiral  slab  for  a  TMz  plane  wave  incident  at  30°  off  the  x-axis. 
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Figure  4.14:  Internal  fields  along  y  =  0  for  the  inhomogeneous  chiral  slab  shown  in 
Figure  4.13. 
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chiral  materials  in  Figure  4.14.  Figure  4.15  shows  the  inhomogeneous  chiral  slab 
of  Figure  4.13  with  a  TE z  incident  plane  wave  from  30°  off  the  x  —  axis.  The  most 
significant  change  from  Figure  4.11  is  that  the  field  scattered  outside  the  main  broad¬ 
side  beam  is  at  least  5  dB  higher  for  both  co-  and  cross-polarized  fields.  Comparing 
Figures  4.13  and  4.15,  the  cross-polarized  scattered  fields  are  significantly  different, 
unlike  in  previous  comparisons  between  the  cross-polarized  fields  for  TMz  and  TEz 
incident  plane  waves.  This  means  that  Equation  (3.34c),  which  showed  that  the 
cross-polarized  bistatic  echo  widths  of  circular  cylinders  for  TMz  and  TEz  incident 
plane  waves  are  identical,  is  only  true  for  circular  cylinders.  However,  Figures  4.16 
and  4.17  show  that  the  cross-polarized  backscattered  fields  are  identical  for  TMz 
and  TEz  incident  plane  waves.  The  numerical  results  used  to  generate  Figures  4.16 
and  4.17  matched  within  0.003  dB  and  0.01°  for  magnitude  and  phase,  respectively. 


4.5  Half-Plane  Numerical  Results 

This  section  presents  numerical  results  for  TMz  and  TEz  scattering  by  two  geome¬ 
tries,  a  chiral  slab  on  a  half- plane  and  a  double- wedge  covering  the  tip  of  a  half- plane. 
All  data  is  at  a  frequency  of  300  MHz.  The  chiral  parameters  used  in  this  section 
obey  the  limit  set  in  [102]  of  |£r|  < 

Figure  4.18  shows  the  backscatter  echo  width  pattern  for  a  TMz  plane  wave 
(El  =  e-iMrcos^+vs.n*'))  incjdent  upon  a  lossless  chiral  slab  at  the  tip  of  a  perfectly 
conducting  half-plane.  The  slab  is  1  meter  wide  and  0.2  meters  thick,  with  parame¬ 
ters  of  er  =  4.0,  fir  =  1.5,  and  (c  =  0.002  A/V.  For  comparison,  the  echo  widths  of 
the  achiral  slab  and  of  the  bare  half-plane  are  also  shown.  Although  the  half-plane 
is  the  dominant  scatterer,  the  presence  of  the  chiral  slab  does  produce  a  significant 
change  to  the  echo  width.  In  particular,  the  chirality  produces  a  cross-polarized 
component  to  the  echo  width  which,  in  the  region  <  120°,  is  of  comparable  mag¬ 
nitude  to  the  co-polarized  echo  width.  This  cross-polarized  field  is  a  direct  result  of 
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Figure  4.15:  The  co-polarized  and  cross-polarized  bistatic  echo  width  of  a  0.1  m  x 
2.0  m  inhomogeneous  chiral  slab  for  a  TMp  plane  wave  incident  at  30°  off  the  x-axis. 
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Figure  4.16:  The  co-polarized  and  cross-polarized  backscatter  echo  width  of  a  0.1  m 
x  2.0  m  homogeneous  chiral  slab  for  a  TM^  plane  wave  incident  at  30°  off  the  x-axis. 
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Figure  4.17:  The  co-polar.  "d  and  cross-polarized  backscatter  echo  width  of  a  0.1  m 
x  2.0  m  homogeneous  chiral  slab  for  a  TEz  plane  wave  incident  at  30°  off  the  x-axis. 
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Figure  4.18:  The  co-polarized  and  cross-polarized  backscatter  (4>  =  4> ')  echo  width 
of  a  lossless  chiral  slab  at  the  tip  of  a  perfectly  conducting  half-plane. 


the  rotation  of  polarization  which  occurs  in  a  chiral  medium.  The  magnitude  of  the 
internal  fields  along  the  center  line  (y  =  0.1  m)  of  the  slab  is  shown  in  Figure  4.19  for 
<f>'  =  60°.  Again,  the  cross-polarized  fields  (Ex,  Ey ,  Hz)  are  a  result  of  the  rotation 
of  field  polarization  in  the  chiral  medium,  and  are  of  comparable  magnitude  to  the 
co-polarized  fields  ( Ez ,  Hx,  Hy). 

Figure  4.20  shows  the  bistatic  scattering  from  a  perfectly  conducting  half-plane 
with  a  lossy  chiral  double- wedge  covering  the  tip.  The  double- wedge  is  2  meters  wide 
and  has  a  maximum  thickness  ot  0.08  meters  at  the  center.  The  upper  graph  shows 
the  co-polarized  and  cross-polarized  bistatic  echo  widths  for  a  TM^  wave  incident 
from  180°  ( E\  =  e~jk"r).  The  bistatic  echo  width  for  an  identical  achiral  (£r  =  0) 
double- wedge  and  for  the  bare  half-plane  is  also  shown  for  comparison.  The  lower 
graph  displays  the  same  data  for  a  TE7  wave  incident  from  0°  (Hi  =  e]k''T).  In 
each  case,  the  chiral  double  wedge  causes  a  significant  modification  to  the  scattering 
from  the  half-plane.  For  example,  for  the  TM^  case,  the  chiral  wedge  reduces  the 
edge  on  backscatter  (</>  =  180°)  echo  width  from  -8  dB/m  to  about  -60  dB/m,  while 
the  achiral  value  is  about  -16  dB/m.  The  chirality  also  produces  a  significant  cross- 
polarized  component  for  both  the  TM7  and  TE^  cases.  Again,  these  cross-polarized 
fields  are  a  direct  result  of  the  rotation  of  polarization  in  chiral  medium. 
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Figure  4.19:  Internal  fields  along  y  =  0.1  meters  for  the  chiral  slab  and  half- plane 
geometry  shown  in  Figure  4.18  with  a  TMz  incident  wave  from  60°. 
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4.6  Summary 


This  chapter  presented  a  hybrid  integral  equation  method  of  moments/Green’s  func¬ 
tion  solution  to  the  problem  of  TM^  and  TE z  scattering  by  an  inhomogeneous  chiral 
cylinder  of  arbitrary  cross  section  in  free  space  or  in  the  presence  of  a  perfectly  con¬ 
ducting  half-plane.  The  volume  equivalence  theorem  for  chiral  media  was  used  to 
formulate  a  pair  of  coupled  vector  integral  equations  for  the  equivalent  electric  and 
magnetic  volume  polarization  currents  representing  the  chiral  cylinder.  For  chiral 
cylinders,  this  pair  of  vector  equations  was  shown  to  be  equivalent  to  six  coupled 
scalar  equations,  which  included  coupling  between  the  usual  TM  and  TE  polariza¬ 
tions.  These  equations  were  solved  by  the  method  of  moments,  and  numerical  data 
was  shown  illustrating  the  convergence  and  accuracy  of  the  method  of  moments 
solution.  As  in  the  previous  chapter,  the  chirality  of  the  cylinders  caused  a  signif¬ 
icant  rotation  in  the  polarization  of  the  scattered  fields.  The  result  was  that  the 
cross- polarized  fields,  in  the  near  and  far  zone,  can  be  as  large  or  larger  than  the  co¬ 
polarized  fields.  In  addition,  the  numerical  results  showed  that  the  cross-polarized 
backscatter  echo  width  of  at  least  one  noncircular  geometry  is  identical  for  TM;/ 
and  TE;/  incident  plane  waves.  An  examination  of  the  internal  fields  of  one  chiral 
cylinder  showed  that  the  number  of  cells  needed  in  the  method  of  moments  solution 
increased  by  a  factor  of  5.  Part  of  this  increase  is  due  to  the  short  wavelength  of 
one  of  the  circular  polarizations  in  chiral  media.  The  remaining  increase  remains 
unexplained  and  may  be  due  to  the  interaction  of  the  two  circular  polarizations  with 
different  wavelengths. 
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Chapter  5 

A  Microstrip  Line  on  a  Chiral 
Substrate 


This  chapter  presents  a  spectral-domain  Galerkin  method  of  moments  solution  for  a 
microstrip  transmission  line  on  a  chiral  substrate.  Most  chiral  guided  wave  structures 
have  bifurcated  modes  [103]— [109],  i.e.,  pairs  of  modes  with  the  same  cutoff  frequency. 
The  chiral  microstrip  line  does  not  have  bifurcated  modes;  thus,  the  dispersion 
curves  are  single  valued.  However,  the  longitudinal  fields  and  currents  each  contain  a 
significant  asymmetric  component  which  is  not  present  for  an  achiral  microstrip  line. 
Similarly,  the  transverse  fields  and  currents  both  contain  a  symmetric  component 
which  is  absent  in  the  achiial  case. 

The  propagation  constant  and  currents  for  the  microstrip  line  are  calculated  us¬ 
ing  the  spectral- domain  approach  [130,  131],  which  is  a  Galerkin  method  of  moments 
solution  where  the  impedance  matrix  and  fields  are  calculated  by  numerically  eval¬ 
uating  the  Fourier  transform  of  a  spectral  Green’s  function.  This  Green’s  function 
accounts  for  the  effects  of  the  chiral  substrate  and  the  perfectly  conducting  ground 
plane.  The  chiral  substrate  causes  the  additional  longitudinal  asymmetric  and  trans¬ 
verse  symmetric  field  components,  requiring  the  expansion  for  the  longitudinal  and 
transverse  currents  to  contain  both  even  and  odd  components.  This  requirement 
is  satisfied  by  expanding  each  current  in  terms  of  Chebyshev  polynomials  weighted 
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Figure  5.1:  Microstrip  line  on  a  grounded  chiral  slab. 

by  the  edge  conditions  [131]— [135] .  A  detailed  discussion  on  edge  conditions  can  be 
found  in  [136,  pp.  385-387]. 

5.1  Microstrip  Field  Expansions 

The  geometry  of  the  problem  and  the  field  expansions  for  each  region  are  presented 
in  this  section.  The  open  microstrip  line  on  a  grounded  chiral  slab  is  shown  in 
Figure  5.1,  where  the  slab  (region  II)  has  constitutive  parameters  (/u,  e,  £c)  and 
thickness  T.  The  right  and  left  wave  numbers,  kn  and  ki,  and  the  wave  impedance, 
T)c,  are  defined  in  Equations  (2.12)  and  (2.11).  The  microstrip  line  is  W  wide, 
infinitely  thin,  and  perfectly  conducting  with  a  current  distribution  of  J(x)e~*k'z. 
Region  I  (y  >  T)  is  free  space,  with  parameters  ( fi0 ,  e0)  and  wave  number  ko  = 

5.1.1  Field  Expansions  in  the  Free  Space  Region 

First,  consider  the  fields  in  region  I.  Since  this  is  a  source-free  achiral  region,  these 
fields  may  be  expanded  as  the  sum  of  a  TMy  field  and  a  TEy  field  [115,  sec.  3-12]. 
These  fields  may  be  constructed  from  the  magnetic  and  electric  vector  potentials 
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given  by 
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where  kv  =  —jyjk*  +  k*  —  satisfying  the  wave  equation  and  the  radiation  con- 
dition  (see  Appendix  E).  The  functions  A  and  F  are  the  Fourier  transforms  of  Ay 
and  Fy  with  the  term  e~^kvV+klZ^  factored  out,  and  are  determined  by  enforcing  the 
boundary  conditions  at  y  =  T.  The  TMr  fields  are  obtained  from  Equation  (5.1a) 
using  H  =  V  x  A 
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The  TE)  fields  are  obtained  from  Equation  (5.1b)  using  E  =  —V  X  F 
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5.1.2  Field  Expansions  in  the  Chiral  Slab 

In  the  chiral  slab,  the  fields  are  expanded  in  terms  of  the  right  and  left  circular 
vector  potentials  developed  in  Section  2.3.  To  satisfy  the  boundary  conditions  of 
zero  tangential  electric  field  at  y  =  0,  a  quasi-TMy  field  can  be  formed  by  the  sum 
of  a  RCy  field  and  a  LCy  field.  These  RCy  and  LCy  fields  are  formed  from  the 
vector  potentials 

oc 

rv,m  =  —  J  cosky. nye^I+k^dkx  (5.4a) 

-  oc 
.  oc 

Um  =  j  Qm  COS  kV'Lye-J^+k"’dkx  (5.4b) 

—  OC 

where  =  —jyjkl  +  k*  -  k2R  and  ky_i  —  —j\jkl  4-  k]  -  k2L.  The  resulting  field  is 
a  TMy  field  if  (c  =  0;  hence,  the  name  quas;-TMy  Similarly  »  quasi-TEv  field  can 
be  formed  from  the  vector  potentials 


Rv,e  =  ^  /  <?E^Binfcv,flye-^r+‘" 


]dkT 


“OC 

oc 


Ly,e  =  ^  j  QE~~^mKLye~^I+k^dkx. 


(5.5a) 

(5.5b) 


The  functions  Qm  and  Qe  are  determined  when  enforcing  the  boundary  conditions 
at  y  =  T. 
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The  quasi-TMy  electric  fields  are  calculated  by  substil 
:tor  potentials  from  Equation  (5.4)  into  Equation  (2.24) 

nr  /  \ 


me  quaw-  i  iUy 
vector  potentials  from  Equs 

f 

—  oc 

—  oc 

«■  ■ 

—  OC 

=  s7M*& 

“OC 

p  _  1  In  k*  +  k> 

£i.,y  -  /  Q*t—r— 


ting  the 


£/.,  =  f 

27T  j 


JQM~kT 

—  OC 

7  ««  &r 

-oc 


■  sin  fcy.fi  y  +  fc2  cos  fcy.R  y  j 

e-jikrT+k.z)^ 

(5.6a) 

cos  fcy.fi  y  e~}i'kzX+kz^dkI 

(5.6b) 

sin  fcy.R  y  -  kx  cos  fcy.R  y^ 

|  e-j(kxx+k.z)dkr 

(5.6c) 

-  sin  fcy,/.  y  -  kz  cos  fcu,y,  y^ 

e~^rT  +  k,z)dki 

(5.6d) 

cos  kv,Lye-’(kzX+kzX)dkx 

(5.6e) 

1  sin  fcyX  y  +  fcr  COS  fcy.L  y^ 

e-H***+k'*)dkr' 

(5.6f) 

-oc 

Similarly,  the  quasi-TEy  electric  fields  are  calculated  by  substituting 
LCj  vector  potentials  from  Equation  (5.5)  into  Equation  (2.24),  pro 

En.x  =  ~  J  Qe  sin  *v,R  y  ~  k  cos  fcy,fi  yj  e~AhzX+k,z)dkx 

—  OC 

r  _  _L  /  rj  k*  ^  r.ir  XI  p~}(kxT+k‘z^dk 
&n,v  —  7r~  I  Va  i  sin  ft.y  fi  y  e  a«r 

Z7T  y  «„  r 

—  OC 

-oc  '  ' 

£l,t  =  ~  J  Qe  sin  ky  Ly  +  fcx  cos  fcy.z.  y  j  e~}{kzX+k,z)dkT 

-oc  ' 

£.  =  — L  7 QEk*  +  k *  sin  kyLy  e~^kzI+kzz)dkx 

2ir  J  fcu.r 


ing  the  RC> 


(t~^  S‘n  ^  V  + 

-oc  '  V’ 

1  7n  kl  +  V  ■  i. 

-  —  /  Qe~t - sinfcuxye 

27T  y  fcyX 

—  oc 
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£/.,  =  ^  j  Qe  (-^y  sinfcI/.Ly  +  k,coskV'Ly\  e-«k'r+k">dkx.  (5.7f) 

—  C5C 

The  magnetic  fields  can  now  be  determined  using  Equation  (2.18)  of  Section  2.2 


(  H «  ^ 

■~»l 

i 

(  En  ' 

\n,  ) 

Vc 

l  ~E'-  , 

5.2  Coefficients  of  the  Field  Expansions 


This  section  presents  the  solution  for  the  coefficients  in  the  field  expansions.  Once 
these  coefficients  are  known  the  fields  may  be  calculated  in  any  region.  The  co¬ 
efficients  are  determined  by  enforcing  the  boundary  conditions  at  the  interface  of 
V  =  T: 


K 

=  El' 

(5.9a) 

El 

=  El' 

(5.9b) 

Hi 

=  Hi1  -  Jz(x)e~jk‘* 

(5.9c) 

Hi 

=  Hi'  +  Jr(x)e~jk" 

(5.9d) 

where  Jz(x)  and  Jz(x)  are  the  components  of  the  current  J(x)  on  the  microstrip 
line.  The  fields  in  region  I  {E[,  E[,  H't,  H'z)  are  the  sum  of  the  appropriate  field 
components  from  Equations  (5.2)  and  (5.3).  The  electric  fields  in  region  II  {E'J , 
E^1)  are  the  sum  of  the  appropriate  field  components  from  Equations  (5.6)  and 
(5.7).  The  magnetic  fields  in  region  II  {Hi' ,  H[l)  are  determined  using  the  same 
field  components  and  Equation  (5.8). 

To  solve  the  boundary  condition  equations  using  the  spectral-domain  approach 
[130,  131],  Equation  (5.9)  is  Fourier  transformed  from  the  space  domain  (x)  into  the 
spectral  domain  (/er),  evaluated  at  y  =  T,  and  after  some  simplification,  results  in 

^ A-kzF  =  QM{kTC+ +k2A-)  +  QE(k,B+ -kTA~)  (5.10a) 

u/c„ 
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*4  I + *4 

Wfo 

=  Qm{KC+  -kTA~)  +  QE(-kIB+  -KA-) 

(5.10b) 

jk,A  +  jkrKF 

Qm  ^ 

Vc 

-kxC~  -  k,A+)  4  4  fc,A+)  -  J, 

Vc 

(5.10c) 

jk,A  +  jk‘k-F 

II 

cj  *Oi 

a 

~kzC~  4  JfcTA+)  +  —  (fc,B-  4  fc,A+)  4  J, 

Vc 

(5.10d) 

v/Hcrc 

A±  = 

=  COS  fcy  rT  4  COS  kyjT 

(5.11a) 

= 

~  sin  ku  nT  i  t —  sin  A:,,  /  T 

ky.H  ' 

(5.11b) 

C*  = 

^ysR  •  i  /ti  i  ^y,L  •  »  rp 

=  — — sin  kyj{l  ±— — sin  kyj  l. 
k/{  k/ 

(5.11c) 

The  qualities  Jx  and  J2  are 

the  Fourier  transforms  of  the  components  of  J(x)  and 

are  defined  by 

OC 

Jx, 

r(*)  =  ~J  l,(K)e-jkrTdkx. 

—  rjc 

(5.12) 

Equation  (5.10)  can  be  rewritten  as 

u>e„ 

=  QmC+  -  Qf.A 

(5.13a) 

-F 

—  QmA~  4  Qe&+ 

(5.13b) 

j  *'  F 
W/iu 

r,  C~  ^  A+  kxJz  -  kzJx 

Vc  Vc  K  4  fcj 

(5.13c) 

jA 

=  A+  s  kz  Jz  4  kxJx 

=  Qmv<  Qev<  *,’  +  *? 

(5.13d) 

where,  for  example,  Equation  (5.13a)  was  obtained  multiplying  Equation  (5.10a) 
by  kx  and  Equation  (5.10b)  by  kz,  adding  these  products,  and  dividing  through  by 
k\  4-  k After  eliminating  A  and  F,  Equation  (5.13)  can  be  written  in  matrix  form 
as 

j^c+  4  -jipA-  4 

J  rie  J  kv  nr 


J  U>/i0  T)c  +  W*l0  Vc 


Qm 

1 

Qe 

"  *2  4  fc? 

fax*) 2 

2  X 


(5.14) 


96 


Equation  (5.14)  can  now  be  solved  for  Qm  and  Qe-  in  order  to  find  the  fields  due 
to  the  *  and  z  components  of  the  current,  Q A/  and  Qe  are  written  as 


Qm  =  Q  a/,x  4-  Qm.z 
Qe  =  Q  E,x  +  Qe.z 


where 

Qm.x 

Qm.z 

Qe.t 

Qe.z 


(j  K  B*  +  ) 

,  (  U'e°  A~ 

~  Kz  J-7—  A 

-~)1 

V  w/io  Vc  J 

\  *y 

Vc  )\ 

-4^  +  dlW(AA-_£:)] 

V  K  Vc )  \  <^mo  vc  )\ 

k,  (i-TC*  +  —)+*.  (}—A-  -  — ) 

V  kv  Vc  J  \  u>Hu  Vc  / 


(5.15a) 

(5.15b) 


(5.16a) 


=  £r\~k,  (ji, B*  +  —  )+k,(j^A--—  )|  (5.16b) 

2A  V  Vc  J  \  *y  Vc  } 


(5.16c) 

(5.16d) 

(5.16e) 


V  '  ([.M  Mo  J  1*1/  M  *2  Mo 

5  =  COS  ky,RT  COS  ky  lT  —  G  + 

1  (  *«  *y,L  ,  *y,fl  kL  \  ,  rr  .  ,  rr 

L  fa 

U  =  sin  kv  rT  cos  kv^T  4-  sin  fcy^T  cos  fey.^T 

«K  *L  r 

sin  *y,RT  cos  *y,£,T  4-  sin  kv^T  cos  *y,/?T 

*y,R  *y,i. 


U  =  sin  ky  rT  cos  kv^T  4-  sin  kv>iT  cos  kv<RT 
kr  *1, 

V  —  sin  *„,rT  cos  kVyiT  4-  ~r~  sin  K-eT  cos  kv,RT 
ky.R  *y,L 

with  *2  =  W-s/M^c  =  (*fl  +  kc)l 2. 

5.3  Electric  Fields  at  the  Interface 


t  +  ki)j 2. 


(5.17a) 

(5.17b) 

(5.17c) 

(5.l7d) 

(5.17e) 


This  section  uses  the  results  of  the  previous  section  to  calculate  the  Ex  and  E2 
fields  at  the  interface,  which  are  used  in  the  method  of  moments  solution  for  the 
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propagation  constant  k2  and  the  current  J(a:)  on  the  microstrip  line.  These  fields  are 
presented  below  in  terms  of  the  even  and  odd  components  of  the  Fourier  transforms 
of  these  fields  generated  by  x  and  z  polarized  traveling  wave  line  sources  at  the 
center  of  the  strip  (x  =  0,  y  =  T) 

&  -  (5-!8a) 
&  =  +  +  (5.18b) 


Et  -  -Et  =  ’-rrG- 

T]C£S 

U  _  jL.2 

Ei\  =  E{\  =  j— 

7JfZ\ 


(5.18c) 

(5.18d) 


(5.18e) 


where  G  ,  S,  U,  and  F  are  defined  in  Equation  (5.17).  For  example,  the  Ex  field 
generated  by  the  surface  current  Jx(a:)  is  given  by 

=  ^  /  [^(M  +  ££(*„)]  X(M  «-«*■'+*.')  ik,  (5.19) 

—  oc 

where  and  Ej*0  are  the  even  and  odd  components  of  the  Fourier  transform  of 
the  Ex  field  due  to  a  z  polarized  traveling  wave  line  source  at  x  =  0  and  y  =  T. 
In  a  conventional  achiral  microstrip  line  Ex*a,  E^0,  £/*,  and  E^e  are  zero,  causing 
Jz(x)  and  JT(x)  to  be  even  and  odd  functions  of  x,  respectively.  However,  this  is  not 
true  for  a  chiral  microstrip  line,  thereby  requiring  a  set  of  even  and  odd  expansion 
functions  for  each  current  component. 

5.4  Method  of  Moments  Solution 

The  basis  functions  and  impedance  matrix  for  a  Galerkin  method  of  moments  so¬ 
lution  are  presented  in  this  section.  The  A  and  z  components  of  the  current  are 
expanded  as 


Jz(x)e 


-jkjz  _ 


(5.20a) 
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A, 


(5.20b) 


n=0 


where  7X  n  and  I2  n  are  the  unknown  coefficients.  The  basis  functions  Jx  „  and  Jz  n 
are  Chebyshev  polynomials  weighted  by  the  edge  conditions  [131]— [135] 

4  Un{2x/W) 


•M*)  =  ^Tn(2x/W)/y/l-(2x/Wy 


(5.21a) 

(5.21b) 


where  Tn(x )  and  Un{x)  are  Chebyshev  polynomials  of  the  first  and  second  kinds, 
respectively.  A  detailed  discussion  of  edge  conditions  is  presented  by  Van  Bladel 
[136,  pp.  385-387].  From  Appendix  F,  the  Fourier  transforms  of  the  basis  functions 
are 


J  (k)  -  2i"^-1-(^/2-) 

*'  ;  kTW/2 

£,„(*,)  =  jnJn(kxW/  2) 


(5.22a) 

(5.22b) 


where  J„(x)  is  an  nth-order  Bessel  function.  The  first  four  basis  functions  for  each 
current  component  are  shown  in  Figure  5.2,  where  the  functions  have  been  normal¬ 
ized  using  x'  =  x/W.  The  primary  advantage  of  these  basis  functions  is  that  they 
have  closed  form  Fourier  transforms  while  matching  the  edge  conditions.  In  addition, 
although  these  basis  functions  are  not  orthogonal,  they  are  sufficiently  orthogonal 
that  the  7X,„  and  /*,„  of  Equation  (5.20)  are  relatively  insensitive  to  the  number  of 
terms  used  in  the  expansion.  This  is  shown  in  Table  5.1  of  Section  5.6. 

The  coefficients  7X,„  and  /*,„  are  found  by  using  the  method  of  moments  to  enforce 
the  boundary  condition  of  zero  tangential  electric  fields  at  the  surface  of  the  perfectly 
conducting  microstrip  line,  i.e., 

EJ/{x)  +  £/'(*)  =  0  (5.23a) 

EJ/{x)  +  EJ/(x)  =  0  (5.23b) 


for  (a:)  <  W/2  and  y  =  T.  Using  Galerkin’s  method  of  moments,  which  uses  weight¬ 
ing  functions  equal  to  the  basis  functions,  Equation  (5.23)  can  be  written  in  block 
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1/vT 


(ix'3  -  3i')/v 

Figure  5.2:  First  four  b 
components. 


matrix  form  a.s 


Zzx  Zzz 


lx 

Iz 


0 

0 


(5.24) 


where 


o 

ii 

■  Iz, a-2]T-  The  xz 

block 

Exz(0,0) 

.  .  .  ZXz( 0,Arz) 

_  Zxz(Nx, 0) 

•••  2I2(A'IiA,) 

11/2 

1 

tq 

*  o. 

N 

"(*)  JX,m{x)  ejk* 

2dx. 

-11/2 

(5.25) 


(5.26) 


The  term  e}k’z  occurs  in  Equation  (5.26)  because  the  basis  and  weighting  functions 
are  of  the  form 


•^i,n(*  ) 


,-jkzz 


(5.27) 


and  the  inner  product  used  for  the  method  of  moments  solution  is  [118,  sec.  12.2.7] 

11/2 

(w(*),g(*))  =  J  w*(i)-g {x)dx.  (5.28) 

-11/2 

Using  Equation  (5.19)  to  calculate  Ex*,  Equation  (5.26)  can  be  written  as 

2«(*,n)  = - J  EJx‘(kx)Jx,n{kx)Jx<m(-kx)dkx.  (5.29a) 

—  OO 

Impedance  elements  in  the  remaining  blocks  are  obtained  in  a  similar  fashion  and 
are  given  by 


em-k,)z  7 

Zxx(m,n)  = 

- J  Exz(kx)Jx<n(kx)  Jx,m(-kx)  dkx 

(5.29b) 

e:(k;-k.)z  . 

Zzx(m,n) 

- 2tt - / 

(5.29c) 

oo 

7  - , ,  -  - 

Zzz(m,n) 

- 2^ —  /  £/‘(W,n(M^(-M<^x. 

—  oo 

(5.29d) 
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The  values  of  kz  for  which  Equation  (5.24)  can  be  solved  are  the  propagation 
constants  of  the  microstrip  line.  From  linear  analysis,  Equation  (5.24)  can  only  be 
solved  if  the  impedance  matrix  is  singular,  i.e.,  the  determinant  is  zero.  Once  the 
propagation  constant  kz  is  found,  the  coefficients  Jxn  and  Izn  can  be  found.  One  of 
the  most  accurate  techniques  for  this  problem  is  singular  value  decomposition  [137, 
sec.  2.9].  Any  matrix,  Z,  can  be  factored  such  that 

Z  =  UXVT  (5.30) 


where 


UrU =  v1 V = 


(5.31) 


and 


£  = 


in, 


Wn 


(5.32) 


If  a  matrix  is  singular,  at  least  one  in,  will  be  zero.  Then,  the  ith  column  of  V  is  a 
vector  in  the  null  space  of  Z  and  is  a  solution  for  the  vector 


h 

In  practice,  the  in,  will  be  not  exactly  zero,  but  it  will  be  numerically  small. 


(5.33) 


5.5  Evaluating  the  Microstrip  Integral  Equations 

This  section  is  concerned  with  the  efficient  numerical  evaluation  of  the  electric  fields 
and  the  elements  in  the  impedance  matrix.  For  example,  EJS‘  and  Zxz  are  given  by 

i  00 

EJ/(x)  =  —J  EJx'{kx)jt{kx)e-’(k'x+k"Ukx  (5.34a) 

—  oc 

Zxz  = - —  J  EJx‘(kx)Jz(kx)Jx(-kx)dkx  (5.34b) 

—  oc 
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where  the  reference  to  the  particular  elements  of  the  impedance  matrix  has  been 
dropped  for  convenience.  As  in  Equation  (5.19),  the  Fourier  transforms  for  the 
electric  fields  and  currents  can  be  decomposed  into  even  and  odd  components: 

EJx'(kx)  =  +  (5.35a) 

MK)  =  J,Ak*)  +  JzAkA  (5.35b) 

£(-*,)  =  jx,e(fcx)  -  jX,o{K)-  (5.35c) 

Equations  (5.35a)  and  (5.35b)  can  be  combined  as 

EJxzAk*)  =  EtAk*)  +  £i,o(M  /w(**)  (5.36a) 

Tjxzak)  =  EiAK)  JzAkA  +  ££,(*,)  JU*,).  (5.36b) 

Then,  Equations  (5.34a)  and  (5.34b)  can  be  written  as 

°C 

£?'*(*)  =  -  f  EJxzAk x)  coskxx  e~’k'*  dkx 

7T  J 

0 

.  oo 

—  —  [  EJXIAkx)  sinfcxx  e~jkzZ  dkx  (5.37a) 

7T  J 

0 

em-k,)z  °? _ 

Zxz  = - /  EJxtAK )  «/x,e(*x)  dkx 

7T  J 

0 

+ - - - ]  EJx,Akx)  JxAk .)  <^x.  (5.37b) 

o 

Equations  for  the  remaining  electric  fields  and  impedance  elements  can  be  written 
by  changing  the  references  to  x  and  z  in  Equations  (5.37a)  and  (5.37b). 

Numerical  Integration 

Efficient  numerical  integration  of  Equations  (5.37a)  and  (5.37b)  requires  careful  ex¬ 
amination  of  the  kernels  of  the  integrals.  First,  consider  the  region  where  kx  is 
Bmall.  Surface  wave  poles  exist  in  the  region  between  &o  and  the  maximum  of  kn, 
ki,  where  A  of  Equation  (5.17a)  goes  to  zero.  These  singularities  can  be  avoided  by 
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Im(kx ) 


Figure  5.3:  Integration  contour  to  avoid  surface  wave  poles. 

using  the  integration  contour  shown  in  Figure  5.3.  Figure  5.4  shows  the  integrand  of 
the  self-impedance  of  the  Maxwellian  (J2,o)  current  distribution  for  kT  =  0  to  1000, 
a  propagation  constant  of  kz  =  349.3  (Ag/A0  =  0.6)  and  several  values  of  (r  from  0.0 
to  0.005  (the  limit  set  in  [102]  is  j£cj  <  =  0.0053).  The  approximate  locations 

of  ku,  k2,  and  2A:2  are  shown  along  the  Re(kx)  axis  for  refeienv.e,  where  fco  =  209.6 
and  k2  =  419.2  =  u/y'/Ze  for  a  frequency  of  10  GHz,  er  =  4,  fir  =  1,  T  =  3mm,  and 
W  =  3mm.  The  limit  of  |£c|  <  \fe/ fi  corresponds  to  k^i  =  k2(\/2  ±  1).  Due  to  the 
wide  variation  of  the  integration  kernel  shown  in  Figure  5.4,  Romberg  integration 
was  chosen  for  the  integration  contour  shown  in  Figure  5.3. 

Next,  consider  the  region  where  kx  is  large.  The  expressions  in  Equation  (5.18) 
have  the  following  asymptotic  behavior: 


Et% 

oc  kx 

(5.38a) 

Et 

oc  l/kx 

(5.38b) 

Et 

=  ejx:0  a  i 

(5.38c) 

ej,% 

t-H 

8 

N 

1 

II 

(5.38d) 

Et 

II 

i " 

8 
i— * 

(5.38e) 

The  remaining  terms  in  the  integrals  of  Equation  (5.37)  are  sine,  cosine,  and  Bessel 
functions.  Since  the  asymptotic  form  of  a  Bessel  function  consists  of  a  cosine  function 
divided  by  the  square  root  of  the  argument,  all  the  integrals  in  Equation  (5.37)  can 
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/m(lntegrand)  Re(  Integrand) 


Integrand  of  Zn 
Frequency  =  1  0  GHz 
ur  =  1  ,  T  =  3  mm,  W  =  3 


M  kx ) 

Figure  5.4:  Integrand  of  the  self-impedance  of  a  Maxwellian  current  distribution 
(Z„(o,o))  f°r  a  normalized  guide  wavelength  of  0.6  (&,=349.3)  and  a  range  of  chirality 
admittances. 
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be  written  for  large  kx  as 

J  -j^p  c os(uJifcr -f  c«)  cos(u;2 kx  +  (3)  dkx  (5.39) 

o  r 

where  P  =  0.5, 1.5, 2.5  for  electric  field  calculations  and  P  =  2,4  for  impedance 
calculations.  For  achiral  impedance  calculations  P  is  always  2,  because  Ex*a  is  zero 
for  achiral  microstrip  lines.  The  closed  form  of  Equation  (5.39)  is  presented  in 
Appendix  H  in  terms  of  sine,  cosine,  and  Fresnel  integrals.  Since  the  expressions 
of  Equation  (5.18)  approach  their  asymptotic  forms  in  Equation  (5.38)  much  more 
rapidly  then  the  Bessel  function  approaches  a  cosine  function  divided  by  the  square 
root  of  the  argument,  the  use  of  these  asymptotic  forms  is  limited  by  the  asymptotic 
form  of  the  Bessel  functions. 

In  the  intermediate  region  between  2max{kn ,  k /,)  and  the  start  of  the  asymptotic 
region,  the  integrals  in  Equation  (5.37)  are  similar  enough  to  the  asymptotic  forms 
that  the  number  of  integration  points  needed  can  be  predicted  from  the  frequency 
of  the  two  sinusoids  in  the  integral,  i.e.,  |u>i  +  u;2|  and  ju^j  —  u/2|.  An  alternative 
approach  for  integrating  the  self-impedance  of  the  Maxwellian  current  distribution 
(Z„(0,a))  is  presented  in  Appendix  G. 

5.6  Numerical  Results 

This  section  presents  numerical  results  demonstrating  the  accuracy  of  the  method 
of  moments  solution,  and  also  some  effects  of  chirality  on  a  microstrip  transmission 
line.  To  simplify  comparisons  all  results  are  for  a  substrate  of  thickness  T  —  3  mm, 
er  of  4,  and  fiT  of  1,  with  a  microstrip  line  of  width  W  =  3  mm.  Since  the  microstrip 
current  can  only  be  found  to  within  a  constant,  all  currents  are  normalized  so  that 
=  1.  When  plotting  currents  or  fields,  the  real  part  is  shown  as  a  solid  line,  and 
the  imaginary  part  as  a  dashed  line.  The  initial  data  will  illustrate  the  accuracy  of 
the  method  of  moments  solution  by  showing  that  it  satisfies  the  boundary  condition 
of  zero  tangential  electric  field  on  the  perfectly  conducting  strip. 
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The  longitudinal  current  on  a  thin  achiral  inicrosirip  iine  is  essentially  the  Max¬ 
wellian  distribution  1/ \J\  —  (2 x/W)2.  The  transverse  current  is  an  odd  function 
satisfying  the  edge  condition  yj  1  —  (2 x/W)2,  similar  to  the  Jx  \  expansion  function 
of  Equation  (5.21a).  Figure  5.5  shows  the  longitudinal  ( Ez )  and  transverse  ( Ex ) 
electric  fields  at  the  interface  (y  =  T)  when  these  two  expansion  functions  are  used 
in  the  solution  for  thin  chiral  and  achiral  microstrip  lines.  Note  that  this  solution 
produces  essentially  zero  tangential  electric  field  on  the  strip  for  the  achiral  case, 
but  not  for  the  chiral  case.  This  indicates  that  the  method  of  moments  solution 
for  the  achiral  case  is  essentially  correct.  However,  the  relatively  poor  field  match 
for  the  chiral  case  indicates  that  it  is  not  correct,  and  thus  more  expansion  func¬ 
tions  are  needed.  The  real  asymmetric  Ez  field  component  shown  in  the  upper  left 
graph  is  due  to  the  chirality  of  the  substrate.  This  field  component  is  unaffected  in 
the  solution  by  even  Jz  and  odd  Jx  expansion  functions,  thereby  requiring  odd  J2 
and/or  even  Jx  expansion  functions  for  a  reasonably  physical  solution.  The  imagi¬ 
nary  symmetric  E2  and  real  asymmetric  Ex  field  components  for  both  the  chiral  and 
achiral  cases  are  virtually  identical,  indicating  that  modest  chirality  produces  only 
a  relatively  small  perturbation  in  the  currents  and  fields. 

A  more  accurate  approximation  for  the  currents  is  shown  in  Figure  5.6  using  ten 
longitudinal  and  ten  transverse  expansion  functions.  The  left-hand  graphs  show  that 
the  fields  nearly  satisfy  the  boundary  condition  of  zero  tangential  electric  field  on 
the  microstrip  line.  The  corresponding  currents  are  shown  in  the  right-hand  graphs. 
The  even  transverse  current  component,  which  occurs  solely  because  of  the.  chirality, 
is  significantly  larger  than  the  odd  transverse  current  component.  The  calculated 
propagation  constant  for  the  chiral  microstrip  line  using  twenty  expansion  functions 
differs  only  slightly  from  that  of  the  achiral  microstrip  line  shown  in  Figure  5.5.  The 
imaginary  part  of  kt,  corresponding  to  attenuation  on  the  line,  is  negative  and  of 
the  order  10~8. 
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(kV/m)  Ez  (V/m) 


Chiral  &  Achiral  Microstrip  Line 
Two  Lowest  Achiral  Expansion  Functions 
£r  =  4,  jUr=  1,  T  =  3  mm,  W  =  3  mm 
Frequency  =  500  MHz 


Figure  5.5:  Et  and  Ex  fields  at  y  —  T  for  identical  chiral  and  achiral  microstrip 
lines,  solved  using  the  two  lowest  order  expansion  functions. 
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The  variation  of  the  propagation  constant  and  the  current  coefficients  with  re¬ 
spect  to  the  number  of  basis  functions  used  is  shown  in  Table  5.1  for  the  geometry 
and  material  parameters  in  Figure  5.6.  Accurate  calculation  of  the  propagation  con¬ 
stant  requires  very  few  expansion  functions  as  seen  in  Table  5.1.  However,  accurate 
calculation  of  the  imaginary  part  of  k2  requires  an  accurate  representation  of  the 
current  on  the  microstrip  line,  which  is  dependent  on  the  width  of  the  microstrip 
line  in  wavelengths.  For  comparison,  Figure  5.7  shows  the  electric  fields  and  currents 
for  the  microstrip  line  examined  in  Figures  5.5  and  5.6  at  5  GHz. 

The  dispersion  curve  shown  in  Figure  5.8  shows  the  normalized  guide  wavelength 
(Ag/A„)  for  the  fundamental  mode  of  a  chiral  microstrip  line,  for  a  range  of  chiral 
parameters.  The  case  £c  =  0  corresponds  to  an  achiral  line.  As  seen  in  Figures  5.5 
and  5.6,  the  chirality  of  the  substrate  significantly  modifies  the  fields  and  currents  of 
the  microstrip  line  even  for  small  chiral  parameters.  However,  as  Figure  5.8  shows, 
the  propagation  constant  is  not  significantly  affected  unless  the  chiral  parameter  is 
a  significant  fraction  of  its  maximum  value  set  in  [102]  of  ( f,mai  =  y/e/V,  which  in 
this  case  is  0.0053. 

Figure  5.9  shows  the  relationship  between  substrate  thickness  and  normalized 
guide  wavelength  for  chiral  and  achiral  microstrip  lines  at  frequencies  of  0.5  GHz 
and  2  Ghz.  The  relationship  between  microstrip  line  width  and  normalized  guide 
wavelength  is  shown  in  Figure  5.10  for  the  same  frequencies. 

5.7  Summary 

This  chapter  has  presented  a  spectral-domain  solution  to  the  problem  of  propagation 
on  a  microstrip  transmission  line  on  a  chiral  substrate.  Circular  vector  potentials 
were  used  to  expand  the  fields  in  the  chiral  substrate.  The  primary  effect  of  the 
chiral  substrate  is  to  generate  asymmetric  longitudinal  and  symmetric  transverse 
fields.  This  effect  could  significantly  alter  the  properties  of  microwave  and  higher 
frequency  devices  constructed  on  a  chiral  substrate. 


110 


C/5 

G 

O 


U 

es 


S 


,0 


a 


§ 

«  10 
S 


o  §3 


C>h 

u 

v  d 
a  •- 

03  to 
V  Irt 

i> 


-a  4$ 
.-*3  d 

*  3 
“  3 

d  cl, 
U  — ! 

(J  ctf 

c 

to  .2 

o  id 

“  6 
s"S 

G  id 


tfi 


3 


V 

M 


-a 
ts 
a ! 


<u 

B 

o 

bo 


v  43 
(3  JC 

id  *» 


d  0 
o  d 


d  .2 


.8  d 


id  o 

tO  03 
S3  03 

P-,  ■*» 

o  d 
i;  to 
a,  g 

g  o 
£  6 


O  o 
.11 


*»  -d 
ed  — 
•C  13 

J  a 

^  « 
•  ■  _d 


i/5  d 
^ 


«5 


H  3 


in 


(kV/m)  Ez  (V/m) 


Chiral  Microstrip  Line 
er  =  4,  fJLr  =  1,  T  =  3  mm,  W  =  3  mm 
£c=.001,  Frequency  =  5  GHz 
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Figure  5.7:  Fields  and  currents  at  y  =  T  for  a  chiral  microstrip  at  5  GHz  solved 
using  ten  Jz  modes  and  ten  Jx  modes. 
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Normalized  Guide  Wavelength 


Chiral  Microstrip  Line 

£r  =  4,  jJLr=  1,  T  =  3  mm,  W  =  3  mm 


Figure  5.8:  Normalized  guide  wavelength  (A9/A0)  versus  frequency  for  the  funda¬ 
mental  mode  of  chiral  and  achiral  microstrip  lines,  for  a  range  of  chiral  parameters. 
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Normalized  Guide  Wavelength 


Chiral  Microstrip  Line 
£r  =  4,  JUr=1,  T  =3  mm 


Width  (mm) 

Figure  5.10:  Normalized  guide  wavelength  (A9/A0)  versus  microstrip  line  width  for 
the  fundamental  mode  of  chiral  and  achiral  microstrip  lines  at  two  frequencies. 
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Chapter  6 
Conclusions 


In  this  report,  the  interaction  of  objects  composed  of  chiral  media  with  electromag¬ 
netic  field"  was  analyzed  using  eigenfunction  and  method  of  moments  techniques. 
The  historical  basis  for  the  constitutive  relationships  used  for  chiral  media  was  dis¬ 
cussed  in  Chapter  1.  These  constitutive  relationships  were  shown  to  have  a  direct 
relationship  to  the  constitutive  relationships  for  optically  active  media,  which  can 
be  derived  using  quantum  mechanics. 

Chapter  2  presented  one  method  for  solving  the  wave  equation  in  chiral  media. 
Key  parts  of  this  method  were  then  used  to  develop  circular  vector  potentials,  which 
produce  circularly  polarized  fields  useful  for  field  expansions  in  chiral  media.  Next, 
the  volume  equivalence  theorem  was  developed  for  chiral  media  using  the  same 
procedure  as  that  used  for  achiral  media,  and  then  was  used  to  develop  image  theory 
for  chiral  media. 

An  eigenfunction  solution  for  scattering  by  a  multilayer  circular  chiral  cylinder 
with  a  center  cylinder  described  by  an  impedance  surface  was  developed  in  Chapter  3. 
The  essential  difference  from  achiral  eigenfunction  solutions  was  the  use  of  circularly 
polarized  vector  wave  functions.  The  formulation  was  such  that  the  TE^  solution 
differed  from  the  TMz  solution  by  only  a  plus/minus  sign  at  one  step.  Previous  work 
required  completely  different  solutions  for  TMz  and  TE z  incident  plane  waves.  An 
examination  of  the  scattered  fields  showed  that  the  cross-polarized  bistatic  echo 
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width  of  circular  cylinders  is  identical  for  TMz  and  TE 2  incident  plane  waves.  The 
most  significant  result  of  this  chapter  was  the  effect  of  chirality  on  the  internal  fields. 
In  one  case,  more  Jz  than  J0  current  was  induced  on  a  perfectly  conducting  cylinder 
coated  by  two  layers  of  chiral  media  for  a  TE?  incident  plane  wave,  even  though 
the  achiral  problem  would  have  induced  no  J:  current  on  the  conducting  center 
cylinder.  An  examination  of  the  relationship  between  chirality  and  scattered  fields 
found  that  even  through  the  achiral  bistatic  scattering  pattern  may  have  no  deep 
nulls,  the  co-polarized  bistatic  pattern  for  a  chiral  cylinder  can  have  deep  nulls  for 
certain  ranges  of  £c.  In  one  case  the  backscattered  field  was  entirely  cross-polarized 
for  a  small  range  of  the  chirality  admittance. 

A  pulse-basis  point-matching  method  of  moments  solution  for  scattering  by  an 
inhomogeneous  object  composed  of  chiral  media  in  free-space  or  in  the  presence  of  a 
perfectly  conducting  half-plane  was  presented  in  Chapter  4.  The  numerical  results 
showed  that  the  cross-polarized  backscatter  echo  width  of  at  least  one  noncircular 
geometry  is  identical  for  TM7  and  TEz  incident  plane  waves.  However,  the  cross- 
polarized  bistatic  echo- width  can  differ,  which  is  in  contrast  to  the  circular  cylinders 
of  the  previous  chapter.  An  examination  of  the  internal  fields  of  one  chiral  cylinder 
showed  that  the  number  of  cells  needed  in  the  method  of  moments  solution  increased 
by  a  factor  of  5.  Part  of  this  increase  is  due  to  the  short  wavelength  of  the  circu¬ 
lar  polarized  fields  in  chiral  media.  The  remaining  increase  remains  unexplained 
and  may  be  due  to  the  interaction  of  the  two  circular  polarizations  with  different 
wavelengths. 

A  spectral-domain  Galerkin  method  of  moments  solution  for  a  microstrip  trans¬ 
mission  line  on  a  chiral  substrate  was  developed  in  Chapter  5.  The  propagation 
constant  of  the  transmission  line  showed  only  a  minor  change  due  to  the  chirality 
of  the  media.  In  contrast,  the  fields  were  significantly  changed.  These  fields  can 
be  divided  into  the  symmetric  and  asymmetric  field  components.  The  symmetric 
field  components  were  only  slightly  affected  by  the  chirality  of  the  media.  However, 
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the  asymmetric  field  components  were  due  in  entirely  to  the  chirality  of  the  media. 
These  asymmetric  field  components  would  significantly  affect  the  coupling  between 
microstrip  lines  on  a  chiral  substrate. 
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Appendix  A 

Scattering  by  a  Chiral  Coated 
Perfectly  Conducting  Cylinder 


This  appendix  presents  an  eigenfunction  solution  for  scattering  by  a  perfectly  con¬ 
ducting  circular  cylinder  coated  by  a  single  homogeneous  layer  of  chiral  media.  The 
development  is  similar  to  the  more  general  case  of  scattering  by  a  multilayer  circular 
chiral  cylinder  presented  in  Chapter  3.  The  geometry  of  the  coated  cylinder  is  shown 
in  Figure  A.l,  where  the  constitutive  parameters  of  the  chiral  coating  are  (/i,e,£r) 
with  an  outer  radius  of  b  meters.  The  right  and  left  wave  numbers,  kR  and  k l,  and 
the  wave  impedance,  tjc,  of  the  coating  are  defined  by  Equations  (2.12)  and  (2.11). 
The  perfectly  conducting  cylinder  has  a  radius  of  a  meters.  The  external  medium  is 


free  space,  with  parameters  (/x0?e0),  wave  number  and  characteristic 

impedance  t/0  =  o- 


External  and  Internal  Field  Expansions 

This  section  presents  the  eigenfunction  expansions  for  the  fields  external  and  inter¬ 
nal  to  the  chiral  coated  cylinder.  The  external  scattered  fields  are  expanded  as  in 
Equation  (3.5) 

ES  =  E  l  "  (A. la) 
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Figure  A.l:  Geometry  of  a  chiral  coated  perfectly  conducting  cylinder. 

Hs  =  J-  £  j-  kMWfi.,)  +  (A. lb) 

Vi)  n=  —  oc 

where  Nn  and  Mn  are  defined  in  Equation  (3.2). 

Using  the  right  and  left  circularly  polarized  vector  wave  functions  [20]  from  Equa¬ 
tion  (3.6),  the  internal  fields  are  expanded  as 

E  =  £  j-  [cXi,  +  ^Ei'.l  +  h„E £,  +  !„eW]  ( A.2a) 

n~  —  oc 

H  =  J  £  j~"  [cnESJi  -  <(„El'i  +  K eS?.’„  -  i„E£]  .  (A.2b) 

TMz  Incident  Field 

This  section  presents  the  solution  for  a  TM z  plane  wave  normally  incident  from 
<f)  =  180°.  As  in  Equation  (3.3),  the  incident  fields  are  expanded  as 


E‘  =  ze-ik°*  =  £  i-"N(')(M 

n=  — oc 

H'  =  -y-e-'1”'  =  -  £ 

V»  yOn=-ac 


(A. 3a) 


(A. 3b) 


The  unknown  coefficients  are  found  by  enforcing  the  boundary  conditions  of  con¬ 
tinuous  tangential  electric  and  magnetic  fields  at  p  =  6  and  zero  tangential  electric 
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fields  at 

p  —  a.  The  scattered  field  coefficients,  a„  and  b„  are  given  by 

' TM  _ 

Bl"(L)A™  4  Bl!KR)B™  4  BV\L)C™  4  Bl?\R)D™ 

(A. 4a) 

ar,  - 

ATM 

tTM 

°n 

2JV< 

nuipobA™’ 

')  +  -  A'<,’2>t/<2,>  -  . 

(A. 4b) 

The  internal  field  coefficients  c„,d„,hn ,  and  ln  are  given  by 

TM  Vvc 

C"  TUopub  A^A/ 

(A. 5a) 

j tm  2JVr  Aj/1 

"  “  xw^A™ 

(A. 5b) 

,,,  2j,r  D™ 

irupub  A£A/ 

(A. 5c) 

tTM  2JVc  C ™ 

"  iru^ftA™ 

(A.5d) 

where 

-  [wl'KL)A™ 

4-  +  IV<2>(I)C™  + 

(A.6) 

and 

A™  = 

^2,2)K.0,(^)  +  "  ^’W*) 

(A. 7a) 

btm  = 

A'i22»K!"(i)  4  -  xW(L) 

(A.7b) 

c'l"  = 

XH  "VP{R)  -  -  A<''2)V„<')(ii) 

(A.Tc) 

D™  = 

A<M,K<2,(£)  -  ^  »*'>(«)  -  a<2"k,<'|(£) 

( A  •  7d ) 

W^(S)  = 

^Zi-l(ksb)H^’(k,b)  -  Z<’>'(M)ff<2,(M) 

Vi) 

(A.7e) 

v;b>)(s)  = 

-  WsW'lM) 

(A.7f) 
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BW(S)  = 

—  Z^(ksb)J^(k(,b)  -  ZM\ksb)Jn{k»b) 

Vo 

(A.7g) 

v(p.<j)  — 

n 

Z^(kLa)Z^’(kRa)  +  ZM'(k,.a)Z^(kRa) 

(A.7h) 

um  = 

w#(M)  -  zir\kLb)z^\kHb). 

(A.7i) 

TE z  Incident  Field 

This  section  presents  the  solution  for  a  TE z  plane  wave  normally  incident  from 
<f>  =  180°.  As  in  Equation  (3.4),  the  incident  fields  are  expanded  as 
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E'  =  yie-jfcoT  = 

■■  E  r"M!,')(/c„) 

n  ~  —  oc 

(A. 8a) 

H’  =  z  —  e~]k{'T 

=  -  £  i- 

(A. 8b) 

Vo 

Von  =  -oc 

The  unknown  coefficients  are  found  by  enforcing  the  boundary  conditions  of  con¬ 
tinuous  tangential  electric  and  magnetic  fields  at  p  =  b  and  zero  tangential  electric 

fields  at  p  =  a.  The  scattered  field  coefficients,  an  and  bn  are  given  by 

„TF.  lTM 

an  =  °n 

(A. 9a) 

,77,-  A<"(L).4™  +  +  A«\L)CJJ  +  A™(R)DlE 

(A. 9b) 

un  - 

The  internal  field  coefficients  cn,d„,hn, 

and  ln  are  given  by 

clE  = 

2 jVc  BrnE 
nvnub  A%E 

(A. 10a) 

dlE  = 

2 jVc  ATnE 
nufiub  AEE 

(A. 10b) 

hlE  = 

2 jVc  DTnE 

TTUIflub  AEE 

(A. 10c) 

llE  = 

2 jVc  ClE 

7 TV  flub  AlE 

(A.IOd) 
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I 

I 

|  where 

A"  =  -  [v;l')(i)J4”;  +  Vi'\R)B™  +  VP(L)CJ,,:  +  V™(R)Dl' }  (4.11) 

(A. 12a) 
(A. 12b) 
(A. 12c) 
(A.12d) 
(A.12e) 


ATnE  =  XMWM{R)  -  2W}^]  -  XMWW(R) 
bte  =  A-p-a)n^i')(i)  -  2ty{)(fl)  -  x<iJ)ny)(r) 

irK/a 

Cl L  =  A-("myi(Ji)  +  -  Xl'^Wl'^R) 

DV  =  X^'HV^(L)  +  -  A-f 'I Wl'KL) 

TTK^a 

AM(S)  =  ^Z^ksb)Jn(k,b)- Z^ksb)J^k, b). 
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Appendix  B 

Evaluating  the  Half-Plane 
Green’s  Function  Integrals 


The  Green’s  function  for  a  z  polarized  line  source  in  the  presence  of  a  perfectly 
conducting  half-plane  involves  the  integrals  [125,  8.3] 


$.(<*,/?)  =  j 


(u2  -f  /32)'~'/2 


(BA) 


for  i  =  1,2.  This  appendix  presents  new  and  efficient  techniques  for  accurately 
evaluating  these  integrals.  For  a  lossless  ambient  medium  a  and  (3  are  real  numbers. 
The  range  of  a  that  needs  to  be  explicitly  considered  can  be  restricted  to  a  >  0 
using  [138] 

S,(a,(3)  =  5,u(/3)-  5,(-a,/3)  (B. 2) 


where  i  =  1,2  and  the  5,()(/3)  are  defined  by 

S,„w  =S,(- 00,0)  (B.3a) 

Sk(P)  =S2(-oo,/?)  =-^e>^2  (y)  +  (y )  .  (B.3b) 

Figure  B.l  shows  the  real  and  imaginary  parts  of  the  integrand  of  Equation  (B.l) 
for  t  =  1  and  /?  =  1.  The  singularities  caused  by  (3  can  be  seen  in  the  lower  graph 
at  u  =  ±j.  The  branch  cut  caused  by  the  square  root  can  be  seen  in  both  graphs 
starting  at  the  singularity  u  —  —j.  Figure  B.2  shows  the  best  regions  for  each 
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Figure  B.l:  The  real  and  imaginary  parts  of  the  integrand  of  the  half-plane  integral 
for  i  —  1  and  0  =  1. 
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Figure  B.2:  Regions  for  each  method  of  evaluating  the  S,  integrals  in  the  half-plane 
Green’s  function  for  a  i  polarized  line  source. 

method  of  evaluating  the  S,  integrals.  These  methods  are  presented  in  the  following 
sections. 


Method  I:  Numerical  Integration  (\a\  >  1.6,  all  (3) 

This  section  presents  a  numerical  integration  technique  for  Si  and  S2  requiring  far 
fewer  points  then  previous  approaches  [129,  138,  139].  Using  the  change  of  variables 
v?  =  x  and  x  =  a2  — jt ,  where  the  second  change  of  variables  corresponds  to  changing 
the  integration  path  for  x  to  the  contour  x  =  a2  to  —j 00  +  a2  to  00,  Equation  (B.l) 


can  be  rewritten  as 


S,(a,/3)  =  -j 


(a2  —  jt  +  (32)'~1/2  y/a2  —  jt 


(BA) 


where  the  portion  of  the  path  at  infinity  does  not  contribute  to  the  integral.  Fig¬ 
ure  B.3  shows  the  resulting  contour  in  the  complex  u  plane.  These  integrals  are  now 
in  a  form  suitable  for  Gauss-Laguerre  Quadrature,  which  can  be  written  as  [140,  sec 
25.4.45] 


J  f(x)e  zdx  =  ■/(*,)  +  Rn 


(B- 5) 
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\ 

\  / 


\  / 
N  / 

\  / 


re  B.3:  Integration  contour  in  the  complex  u  plane. 
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where  w,  and  x,  are  given  in  [140,  table  25.9]  for  N  =  1, ... ,  10, 12, 15  and  R,\  is  the 
remainder.  In  our  case  a  very  conservative  upper  limit  on  R ^  is  given  by 


Itfwl  <  < 


2a2+<* 


for  Sj 


for  S2 


(B.  6) 


2q4+4N 


Empirically  the  N  required  for  a  relative  error  of  8  in  either  5j  or  S2  is  approximately 
given  by 


N  %  -.1  -  log1u(S(.16  -  .031oglo6) 
|a| 


—  O  I  u  \ -  - o  1  u  / 

■5  -  ^  + 11'6  +  R  +  ,og"' s  (-35 -  r)  }. 


(B.7) 


for  0  <  /32  <  2  and  by 


N  *  A  '  R  [L2  ~  R  +  log"' ' 5  {2  5  ~  R  +  logl" 4  (J8  “  r)  } 


(B.  8) 


for  j32  >  2.  These  equations  were  developed  for  N  —  1, . . . ,  10, 12, 15  and  0.00001  < 

8  <  0.01. 


Method  II:  Exact  Integration  (/3  —  0) 

This  section  presents  the  special  case  of  0  =  0,  where  the  5,’s  can  be  exactly  inte¬ 
grated 

S,(|a|,0)  =  ^25, R)  (B.9.) 

(B.9b) 

where  E\(z)  is  the  exponential  integral  given  by 

7  e-JU 

Ei(Z)=f- — du.  {B.  10) 

J  u 

T 

as  defined  in  [140,  sec  5.1]. 
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Method  III:  Small  Argument  Approximations 

(  \a\  <  1.6, /32  <  7  ) 

This  section  presents  efficient  small  argument  approximations  for  the  integrals  S ) 
and  S2  when  a  and  (3  are  small.  In  this  region  the  integrals  are  evaluated  using 

S,{a,/3)  =  ^5l0(/3)  -  sign  (a)  AS,(0,  |aj,/3)  (5-11) 

where  t  =  1,2,  the  S,u  are  defined  in  Equation  (B.3),  and  the  AS,  are  defined  by 


a 

AS,(0,a,/3)  =  I 


l  («2+/?2r,/2 


du. 


(B- 12) 


The  AS,  can  be  evaluated  using  the  Taylor  series  expansion  for  e  ^  and  inte¬ 
gration  by  parts.  For  ASi,  where  /3  0,  this  gives 


AS,(0,  \a\,/3)  =  a,  In 


a 


+ 


\ 


(|)!  +  i  I  -ih yv^  +  02  (*•«) 


where  ax  and  fc2  are  computed  using  reverse  recursion  and  the  relations 

,2n  -  1 


a„  =  1+j- 


2n2 


(3l  an 


+i 


bn  —  J 


,n  —  1  2 


a  b. 


n 


n  +  1  • 


(B.14a) 

(B.14b) 


The  reverse  recursion  is  started  at  n  =  IV,  where  a/v  =  b n  =  1.  Then  the  relative 
error,  6,  is  approximately  given  by 

4 


6  < 


N  N\ 
(2  N)l 


a 


2N-2 


c  ^  r  02 N 

4N P  ' 

For  AS2,  where  (3  ^  0,  the  same  technique  produces 


(B.15a) 

(B.15b) 


A52(0,  \a\,(3)  =  JC]  In 


a 


+ 


1 


(i)  + 1  r  ,/a^b5  {jd' +  j)  (B16) 
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(B.17a) 


where  c\  and  d\  are  computed  using  reverse  recursion  and  the  relations 


.  2n  +  1  „2 


Cn  "  1  +  J2n(n  +  1) 

.  1  2  J 

di  =  ci  +  j-  a  d2 


0*  c„ 


+i 


n  —  1 


c„  -  j- 


a2  dn+]  for  n  >  1. 


(B.17b) 

(B.17c) 


n(n  +  1) 

The  reverse  recursion  is  started  at  n  =  N,  where  ck  =  d/v  =  1.  Then  the  relative 
error,  8 ,  is  approximately  given  by 


8  < 
8  < 


7V(7V+  1)! 

(2  N  +  1)!  2A- 

4 KN\\N  +  1)  P 


(B.18a) 

(B.18b) 


Due  to  the  relative  size  of  the  terms  in  Equations  (B.13)  and  (B.16),  Equa¬ 
tions  (B.15a)  and  (B.18a)  provide  very  accurate  estimates  on  8  for  a  given  a;  how¬ 
ever,  Equations  (B.15b)  and  (B.18b)  give  very  conservative  estimates  on  8  for  a  given 
/3.  Empirically,  the  N  required  for  a  relative  error  of  8  in  S,  is  approximately  given 

by 


N* 


0.85  -  12.3  50  N  +  7.70  a  for  5, 
2.0  -  12.0  8Q  U  +  7.41  a  for  S2 


where  Nm  is  rounded  up  to  N  =  4,8,  or  13  (i.e.  if  N‘  —  4.1,  then  N  =  8). 


(*■19) 


Method  IV:  Small/Large  Argument  Approxima 
tions 


(  H  <  1.6,  P2  >  7  ) 


This  section  presents  efficient  small/large  argument  approximations  for  the  integrals 
Si  and  S2  when  a  <  0  and  0  is  large.  Again  Equation  (B.ll)  is  used  as  in  the  previ¬ 
ous  section.  Then  AS,  is  evaluated  using  the  binomial  expansion  for  (u2  4-  02)^2~' 
and  integration  by  parts.  For  ASi  this  gives 


ASi(0,q,/3) 


1 

\0\ 


eiF(a) 


j_a_ 

4  02 


{B.  20) 
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(B.21a) 


where  e\  and  /2  are  computed  using  reverse  recursion  and  the  relations 

,  .  (2n  —  l)2  1 

e"  ”  1+J  4t>  02  e"+1 

2n  —  1  a2 

fn  =  6n  2n  T2  fn+l 


(B.21b) 


and  F(a)  is  the  Fresnel  integral  given  by 


a 

F(a)  =  j 


e-]U  du. 


{B.  22) 


The  reverse  recursion  is  started  at  n  =  N ,  where  e/v  =  /a’  =  1.  Then  the  relative 
error,  6,  is  approximately  given  by 

(2JV)!2  1 


6  < 


42/v  N\3  02N 


4  (2iV)!  a2A'~2 
<  3  N\2  02"-2' 

For  A S2  the  same  technique  produces 

1 


AS2(0,a,/3)  = 


I0I3 


x  ,3ft.  _ 

9\  F{a)-j~  —  h2  e 


where  g\  and  h2  are  computed  using  reverse  recursion  and  the  relations 

,  .(2n- l)(2n  +  l)  1 

gn  -  l+J  An  p2  9n  +  i 

2n  +  1  a2 

hr,  -  gn  -  —2—  -jp  ftn  +  1. 


(B.23a) 

(B.23b) 

{B.  24) 

(B.25a) 

(B.25b) 


The  reverse  recursion  is  started  at  n  =  N,  where  <7a ■  =  h^  =  1.  Then  the  relative 
error,  6 ,  is  approximately  given  by 

(2N  4-  1  )(2  AT)!2  1 


S  < 


A2NNP 


02N 


4  (2iV  4-  1)!  a2A_2 


(B.26a) 

(B.26b) 


5  4A’iV!2  /?2A'-2 

Empirically,  the  N  required  for  a  relative  error  of  6  in  5,  is  approximately  given  by 

3.2  -  20.8  60-225  +  23.5  $  for  5, 


N‘ 


6.1  -  24.4  60136  -f  20.0  j;  for  S2 


{B.  27) 
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where  TV*  is  rounded  up  to  TV  =  4,8,  or  13  (i.e.  if  TV*  =  4.1,  then  TV  =  8). 

In  summary,  this  appendix  presented  new  and  efficient  techniques  for  accurately 
evaluating  the  integrals  needed  for  the  half-plane  Green’s  function  of  a  line  source. 
For  a  given  relative  accuracy,  these  approximations  are  much  faster  then  the  tech¬ 
niques  presented  in  previous  work  [129,  138,  139). 


132 


Appendix  C 

Chiral  Slab  and  Normal  Incident 
Plane  Wave 


This  appendix  presents  the  solution  for  a  plane  wave  normally  incident  on  a  chiral 
slab  in  free  space,  and  in  a  form  simpler  than  available  in  [54]  and  [62].  Consider 
a  right  circularly  polarized  wave  normally  incident  on  the  chiral  slab  as  shown  in 
Figure  C.l.  The  incident,  reflected,  internal  and  transmitted  fields  are  given  by 


E'  =  {x-jz)ejkny  (C.la) 

Er  =  R{x-jz)e~jkoy  (C.lb) 

E'1  =  T,nt  (*-  jz)eik*y  (C.lc) 

Er‘  =  R,nl(x-jz)e-jk (C.ld) 

Ef  =  T(x-jz)  e’hy  (C.le) 


where  k u  =  Vy/flot ii,  k2  =  &/?,l  =  K  ±  fcj,  kc  =  k ^  —  u;/i£c,  and 

ec  =  e  +  fi£r2-  Simple  matching  of  fields  at  the  boundaries  produces 


Tint 

Rint 


T 

R 


rp  V2  +yc 
2V2 

] i  V2  ~~  VC 
2772 

t)2  (  ejk°'  +  R  e~ik°l 
e*ktl  \  772  cos  kcl  +  jrjc  sin  kcl 
Veff  ~  Vo  ej2k0l 
Vcff  +  Vo 


(C.2a) 

(C.2b) 

(0.2c) 

(C.2d) 
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Figure  C.l:  Geometry  for  an  infinite  chiral  slab 
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Appendix  D 

Fields  of  Circular  Vector 
Potentials 


This  appendix  presents  the  fields  produced  by  the  circular  vector  potentials  devel¬ 
oped  in  Section  2.3.  For  the  *  right  circular  vector  potential 

R  =  sc  (asin&y.H  y  4-  6cos  y)  (D.l) 
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Ey  =  ~j 


( ak2  +  b^xy~ yX  )  sin  kyL  y  +  (bkz  -  a— ~^)  cos  ky>L  y 


,-j(kxX+kzz) 


Ez  = 


(a^--  +  sin  kyj.  y  +  (b^~^  -  akytL)  cos  fcyx  y 

*z.  «/, 


-j(kxX+kzz) 


For  the  y  right  circular  vector  potential 

R  =  y  (a  sin  fcyR  y  +  fc  cos  ky,n  y)  e_',(fclX+*r2) 


(D.4b) 


(D.4c) 


(P.5) 


the  field  components  are 
Ex  =  j 


(akz  +  b^Tyy  R  )  sin  kyM  y  +  (bkz  -  J^~)  cos  fcy,K  V 
kr  kr 


-j(kzx+kzz) 


k2  +  k2 

E„  =  - — -  (a  sin  ky<R  y  +  b  cos  ky  r  y)  e~}^kxT+k’^ 


Ez  =  -j 


( akx  -  b—~kz)  sin  fcy,R  y  +  (fefcT  4-  cos  y 

««  K/i 


,~j(kzx+kzz) 


For  the  y  left  circular  vector  potential 


L  =  y  (a  sin  ky<L  y  +  b  cos  kyj .  y )  e  J(fc*x+**2) 


(D.6a) 

(D.6b) 

(D.6c) 

(P.7) 


the  field  components  are 

Ex  =  j 


( akz  -  b^T^y  L  )  sin  ky,L  y  +  ( bkz  +  a ^  ' )  cos  ky_L  y 

kL  ki 


-}{kzx+kzz) 


L2  j.  LJ 

Eu  = - - -  (a  sin  ky,L  y  4-  5  cos  fcy.L  y)  e“J('CrX+^2) 


P,  =  -J 


ky.L .  L  ..  ,  /tL  _ky,L 


(akx  4-  fc-p^A^sin  fcyX  y  4-  (feZe*  -  a—^kz)  cos  fcy,/,  y 

ki 


-j(kTx+kzz) 


For  the  z  right  circular  vector  potential 

R  =  z  (a sin  &y  «  y  4-  6cos  kyjt  y)  e-J(kiT+k*z'> 


(D.8a) 

(D.8b) 

(D.8c) 

(C.9) 
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the  field  components  are 

k  k  k  k 

Ex  =  -  (a-f— 1 -  +  bky^n)  sin  kyj{  y  4-  (b ■-*—  -  akVtR)  cos  kyHy 

kn  kfi 

e-j(klX+klz)  (D-10a) 

Ey  =  j  ( akx  +  b^~- )  sin  kyJ{  y  +  (bkx  -  a  ky^-z- )  cos  ky  R  y 
l  kn  kr  ' 

e-j(kTx+kzz)  (D,10b) 

L2  i  k2 

Ez  -  (asinky'f{y  +  bcoskytRy)e-j{~hlT+k**\  (D.lOc) 

kR 

For  the  z  left  circular  vector  potential 

L  —  z  (a  sin  ky<i  y  +  b  cos  kytR  y)  e~j(k*x+k‘z)  (£M1) 

the  field  components  are 

k  k  k  k 

Ex  -  (a— p- ^  -  bky  L)sm  ky  L  y  +  (h-f-1  +  aky_L)  cos  ky,L  y 

*L  kL 

e~j(krT+k.z)  (D.12a) 

Ey  =  j  {akT  ~  bt~^-)smky,Ly  +  (bkx  +  -)  cos  kViLy 

ki  ki_ 

e-j(kzx+k,z)  (D12b) 

L2  i  L2 

Ez  = - — — —  (a  sin  ky,L  y  +  bcosky^L  y)  e~j('krX+k‘*\  (D.12c) 

ki 
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Appendix  E 

Spectral  Domain  Fields  in  a  Half 
Space 


This  appendix  examines  the  wave  numbers  of  a  spectral  domain  field  in  a  half  space 
of  y  >  0.  Assume  a  spectral  domain  field  given  by 


OC 


y>  =  J  1 j>{kT)e-Hk'x+k»y+k‘z)dkT 

-oc 

(E.  1) 

where 

fco  =  *£  +  *;  +  *?. 

(E.2) 

For  4-2  propagation,  with  k:  =  f32  —  jat. 

e-,kr'  = 

{E.2) 

which  requires  that  both  a,  and  f3z  be  greater  than  zero. 

y  — »  oo,  therefore  ay  >  0,  using  ky  —  0y  —  jay  and 

The  fields  must  decay  as 

e~JkvV  —  e~avVe~}avV  % 

(EA) 

Take  fcn  and  kT  to  be  real,  then  from  Equation  (E.2) 

K  =  kl  +  01-al  +  fi-al 

(E.5a) 

0  =  ~j2ay0y  -  j2a.0.. 

(E.5b) 
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From  Equation  (E.5b) 


therefore 

0y<O 

and 

-  Jkl  -  kl-k]  for  k2x  +  kj  <  k2 
-j\Jk2  +  k]  -  kl  for  k\  +  k]  >  k2. 


(E.6) 

(E.  7) 

(E- 8) 
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Appendix  F 

Fourier  Transforms  for  Microstrip 
Problems 


This  appendix  presents  some  of  the  Fourier  transforms  used  in  microstrip  problems 
and  how  they  are  obtained.  The  Fourier  transform 


OC 

I  J„( a*)  e-jbx  dx 

-OC 


0  for  |h|  >  a 


(F.l) 


where  T„  is  a  Chebyshev  polynomial  of  the  first  kind,  can  be  obtained  directly  from 
[141,  sec.  6.671.8-10]. 

The  Chebyshev  polynomials  of  the  first  and  second  kinds  have  the  following 
properties 


dd 


Tn(cos6)  =  cos  nd 
sinn# 

£/„_,(  ccs0)  =  — — 

sin  8 

[sin^C/n-^cos^)]  =  nT„(cos0). 


As  a  consequence  of  Equation  (F.2c) 


dx 


VT 


Starting  with 


J„{at)  for  \w\  <  a 

V  a*  —  ur 


(F.2a) 

(F.2b) 

(F.2c) 


(F.3) 


(F.4) 
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use  Equation  (F.2c)  to  get 


_ ey  1 

J„(at )  <-—>  (— j)n —  —  [v/a2  —  w2  U„-\ (u>/a)|  for  |u;|  <  a.  (F.5) 
na  aw  L  J 


The  derivative  property  for  Fourier  transforms  is 

(-Mm  ~ 


(F.6) 


Comparing  Equation  (F.5)  and  (F.6)  the  following  relationship  can  be  obtained 

l  ^[yf^SUn^(w/aj\  for  |u>|  <  a.  (F.7) 

t  na  aw  1  3 


which  can  be  written  as 


J„(ax ) 


J  Ju[o 


e~jbxdx  = 


{-jY-'±VJr=VUn^{bla)  for  |fc|  <  a 
0.  for  |6|  >  a 


(F.8) 


Edge  Modes 


Edge  modes  on  microstrip  lines  can  be  represented  by 


'a  w 

fa  4-  w 


'a2  —  w2  Va2  —  w2 
1  w 

/a2  —  w2  y/  a2  —  w2 


(F.9a) 


(F.9b) 


From  Equation  (F.4),  the  Fourier  transforms  of  these  edge  mode  functions  are  given 


-  [Ju(a<)  -  jaj\(at)\ 


for  It^l  <  a 

Va+1»:  I  I 

0  for  |u>|  >  a 


1  ;  for  |  |  <  a 

-  \  Ju{at)  +  jaJi(at))  < — >  <  Va 


0  for  M  >  a. 


(F.lOa) 


(F.lOb) 


Fourier  transforms  for  these  and  more  complicated  edge  mode  type  functions  can 
also  be  calculated  using  Equations  3.384.1  of  [141]. 
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Appendix  G 

Numerical  Integration  for  the 
Self-Impedance  of  the  Dominant 
Microstrip  Current 


This  appendix  presents  a  technique  which  can  be  applied  to  the  integration  of  the 
self-impedance  of  the  Maxwellian  current  distribution  on  a  microstrip  line.  This  is 
the  term  of  Section  5.4  and  can  be  written  as 


where  f(kT)  approaches  a  constant  for  large  kT. 

Consider  the  integral 

OC 

/  =  J  f(x)dx 

o 

where  f(x )  o c  x~a  for  a  >  1  and  large  x.  The  integral  can  be  written  as 


(G- 1) 


(G- 2) 


/  =  /„(*)  + /,(*)  +  /«(*) 


(G- 3) 


where 


v 

Io(X)  =  j  f(x)dx  (G.4a) 

u 

A’+A 

It(X)  =  J  f(x)dx  (G.4b) 

Y 
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oc 

/«,(*)  =  J  f(z)dx. 


(G.4c) 


Next,  define  a  constant  c  such  that 


-A 

Ii{X)  =  c  J  x- 


(G.b) 


For  large  a: 


_ wn _ 

(a-l)[X»-“-(X  + A)'-]‘ 


Ioo(X) 


oc 

/ •-< 


(G-.6) 


(G.7a) 


*  _ *<■*>,  .  (G.7b) 

(1  +  A/A')  1  -  1  1  ’ 

This  technique  is  used  by  integrating  to  a  point  Xu,  integrating  the  additional 
distance  A,  and  using  Equation  (G.7b)  to  estimate  /<*,( A’o)  and  repeating  for  X\  — 
Xu  +  A.  By  comparing  Ioc{X u)  and  /^(Xi )  an  estimate  of  the  error  can  be  obtained. 

For  the  microstrip  impedance  problem,  A  would  be  chosen  as  the  length  of  one 
or  more  cycles  of  the  sinusoid,  in  this  case  A  =  n  x  and  the  first  X  would  be 
chosen  to  be  near  a  null.  Foi  ,his  special  case  of  a  =  2,  /^(X)  is  given  by 


70O(X)  =  /,(X): 


(G.  8) 
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Appendix  H 

Generalized  Sine  and  Cosine 
Integrals 


This  appendix  presents  generalized  sine 

the  sine  and  cosine  integrals 

and  cosine  integrals  which 

are  extensions  of 

5i(z)  = 

f  sin  t  , 

It* 

0 

(H.la) 

Ci{z)  = 

cos  t 

-/  ,  dt 

2 

(H.lb) 

si{z)  = 

OC  . 

/■  sin  t  , 

- It dt 

2 

(H.lc) 

- 

*  (*)  ~  \ 

(H.ld) 

and  the  Fresnel  integrals 

z 

S(z)  =  /sin  — t 2  dt  =  ~=  /  da 

y  2  y 

o  v  o  v 

-z'1 

z  1  J 
.  /  7T  ,  .  1  f  COS  S  , 

C(z)  =  /  cos  — t2  dt  =  ■■  . —  /  — 7=-  da 

v  '  y  2  y 

0  0 

as  defined  in  [140,  sec.  5.2,  7.3]. 


(H.2a) 

(H.2b) 
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The  first  extension  is  for  higher  powers  in  the  denominator 


f  sint  ,  1  [smz 

si(z,p)  =  -  j  dt  = -pzr^iTr —  +  -  i)J 


OC 

Ci(z,P)  =  - 1 


cos  t  ,  1 

—  dt  = 


p  -  1  It1’ 


cosz  ■ I  P  n 

JpZT  ~  ai(z,P-  1) 


where  si(z ,  1)  =  si(z )  and  Ci(z,  1)  =  Ci(z). 

The  second  extension  is  to  double  cosine  and  sine  functions 

00  1 

CCi{z,P,u\,ot,LJ 2,(3)  =  -  f  jj;  cos(u {t  +  a)  cos(u2t  +  0)  dt 

Z 

^  1 

CSi(z,P,Vi,a, u>2,0)  =  -  J  jp  cosluyt  +  a)sin(u>2t  +  0)  dt. 


2 


Each  of  these  functions  can  be  written  as  the  sum  of  four  terms 

CCi(z,  P, u*i ,  a, a>2, 0)  =  CCi\  +  CCi2  +  CC13  +  CCij 
CSi{z,P,Ui,a,u2,0)  =  CSU+CSi2  +  CSi*  +  CSi. , 


(H.3a) 


(H.3b) 


(H.4a) 


(H.4b) 


(H.5a) 

(H.5b) 


where 


CCi]  — 

COs(a2  ^|w,  w2|f_1  Ci(\u>i  0*1*,  P) 

(H.6a) 

CCi2  — 

c°s(a  +  /3)  ^  +a,!|P-i  Citlw,  +  <*|*,P) 

(H.6b) 

CCi 3  — 

- Sin— ^ — —  |u>j  —  u;2 1 r_  1  -u;2|z,P) 

(H.6c) 

CCiA  = 

~sin(a^+  0)  +a,2| r-i  +  u/2| z,P) 

(H.6d) 

CSn  = 

_cos(a^ — 0)  ^  _a,2|f-i  s^|u,1  -a;2| Z,P) 

(H.6e) 

CSi2  = 

cos(a  +  0)  +U,2|P-J  +  w2|z,  P) 

(H.6f) 

CSi3  = 

_sm(a_0)  |u;i  _W2| r-,  CiQ^  -  u,2\z,  P) 

(H.6g) 
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Sm(y  0)  |w,  +a;2|p-'  Ci(\u>i+u>2\z,P) 


(H.6h) 


CSU  = 


with  adjustments  for  special  cases: 

u>\  —  u>2  <  0:  change  sign  on  CCij  and  CSi\ 
u>i  +  u>2  <  0:  change  sign  on  CCij  and  CSi2 
u)\  —  u>2  —  0:  CCij  -  -  0,  C Six  —  0,  and 


CCi\ 

CSi.x 


cos(o  —  P)  21  p 
2  P  -  1 

sin(a  —  P)  z,_p 
2  P  -  1 


a;i  +  u>2  =  0:  CCiA  =  0,  CSi2  —  0,  and 


CCi2 

CSi* 


cos(a  4-  p)  z'  r 
2  P  -  1 

sin(a  +  /?)  z1_p 

2  p  -  r 


(H.7a) 

(H.7b) 

(H.8a) 

(H.8b) 
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Appendix  I 

Evaluation  of  the  Microstrip 
Impedance  Elements  in  the 
Spatial  Domain 


This  appendix  presents  spatial  domain  equations  and  sample  numerical  results  used 
to  verify  elements  of  the  method  of  moments  impedance  matrix  presented  in  Chap¬ 
ter  5.  Since  the  spatial-domain  Green’s  function  for  a  microstrip  line  does  not  exist 
in  closed  form,  a  substrate  with  material  parameters  (jiu,e0) ls  necessary  to  evaluate 
the  impedance  matrix  in  the  spatial  domain  for  simple  tests.  The  resulting  problem 
can  be  decomposed  into  the  self-impedance  of  a  strip  in  free  space,  and  the  mutual 
impedance  of  the  strip  and  its  image.  The  geometry  of  the  equivalent  problem  is 
shown  in  Figure  1.1. 

The  first  test  involves  numerical  integration  of  the  self-impedance  of  Jz  current 
modes.  The  self-impedance  of  a  Jz  current  mode  in  free  space  is  given  by 

o  "72  1172 

2  __  L2  '  ' 

n 


Zzz,S  — 


4u>e 


^  j  I  H^2\kt\x  -  x'\)  J2{x')dx'  Wz{x)dx  (7.1) 


-ir/2  -ir/2 

and  the  mutual  impedance  of  the  current  mode  and  its  image  is  given  by 
u/2  ir/2 


Zzz,i  = 


4u>e 


z  [  J  HP(k,y/(z  -  x-y  +  it3)  J,(x')  ix'  W,(x)  dx 


(1.2) 


-tr/2  -u/2 

using  Equation  (J.15b)  and  the  Zzz  equivalent  of  Equation  (5.26). 
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Figure  1.1:  Equivalent  geometry  for  self-impedance  of  Jz  microstrip  current  mode. 


For  W  =  0.1,  T  =  0.1,  k0  =  10.479  (500  MHz),  kz  =  10fc„,  and  Jz{x')  = 
Wz(x)  =  1/W,  numerical  integration  of  Equation  (1.1)  plus  Equation  (1.2)  using 
Gaussian  Quadrature  with  280  and  288  points’  results  in  a  value  of  Zzz  —  — j  17663, 
which  differs  from  the  spectral  domain  value  of  Zzz  =  — j  17596  by  0.4  percent.  For 
Jz{x)  =  Wz{x)  —  l/(7r^/W2/4  —  x2),  the  same  numerical  integration  results  in  a 
value  of  Zzz  =  —j  18849,  which  differs  from  the  spectral  domain  value  of  Zzz  = 
—j  19012  by  0.9  percent. 

The  second  test  involves  closed  form  integration  of  the  self-impedance  of  Jx 
current  modes  using  the  small  argument  approximations  for  Hankel  functions.  The 
i-component  of  the  electric  field  at  the  surface  of  the  strip  shown  in  Figure  1.1  due 
to  a  Jx  surface  current  is  given  by 


ir/2 

EX(JX)  =  ~  e-jk“  [  k2z  H™(klP)  +  -  H\2\ktP)  JT(x')  dx' 

4u;e  J  p 

1 1  • /o  r 


using  Equation  (J.13a),  y  =  y',  and  letting  T  ■-*  oo  so  that  the  image  current  does 

not  contribute  to  the  electric  field  at  the  strip.  For  kz  —  0  and  Jx(x')  =  1  /W,  then 

'The  logarithmic  singularity  of  the  Hankel  function  is  avoided  by  using  a  different  number  of 
points  for  the  each  of  the  double  integrals. 


Equation  (1.3)  reduces  to 


-  4=^  7  irbi  ""’w*  -  ■’»  (/-4) 

—  U/2 

Using  the  small  argument  approximation  for  the  Hankel  function  and  integrating, 
Equation  (1.4)  is  approximately  given  by 

For  point-matching  with  kz  =  0,  the  self-impedance  of  a  Jx  current  mode  is 
approximately  given  by 


E  (j )  ~  _  i  .  _ 1 _ L__ 

A  x)  ~  8  27rwelU  [x  -  VF/2  x  +  fV/2 


—  —Et(Jt)  x_0 


J’2 

7ru>eVF2 


For  a  frequency  of  500  MHz,  VF  =  0.06,  and  T  — >  oo,  Equation  (1.6)  gives  a  result 
of  ZTT'S  =  493.5  —  j6358,  compared  to  Zxx,.«  =  465  —  j‘5449  using  spectral  domain 
integration. 

For  a  weighting  function  in  the  shape  of  a  triangle  given  by 


For  a  frequency  of  500  MHz,  W  =  0.06,  and  T  — *  oo,  Equation  (1.8)  gives  a  result 
of  Zxx,»  =  493.5  -  j'8813,  compared  to  Zxx,s  =  490  —  j 7976  using  spectral  domain 
integration.  For  T  =  0.1,  the  mutual  impedance  of  the  image  can  be  calculated  by 
numerically  double  integrating  Equation  (J.13a)  multiplied  by  the  weighting  func¬ 
tion,  resulting  in  Zxxx  =  101.4  -I-  j‘528.2,  compared  to  Zxx,,  =  102  +  j 529  using 
spectral  domain  integration  by  subtracting  the  impedances  for  T  — ♦  oo  and  T  =  0.1. 
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For  k2  ^  0  and  J,(x')  =  1  /W ,  Equation  (1.3)  reduces  to 

M/2 

kl  I  H[?\kt\x  -  x'\)dx' 


EX(JX)  = 


-1 


4weVF 


L  -n/2 


U/2 

~  *'l) d* 


—  U/2 


(1.9) 


Using  the  small  argument  approximations  for  the  Hankel  functions  and  integrating 
Equation  (1.9)  with  the  triangular  weighting  function  given  in  Equation  (1.7),  the 
self-impedance  is  approximately  given  by 


1 


AiUjjeW 


k]  I  W-j-W 


-ir  +  5b2  +  1"Tw’ 

W 


) 


/  w  16  \l 

+fc'(*'T  -Wn2)l  ,U0) 


where  7  =  1.781. 

For  a  frequency  of  500  MHz,  W  =  0.06,  T  — >  00,  and  k2  =  0.1fc();  Equation  (1. 10) 
gives  a  result  of  ZXX3  =  498.4  —  j8817,  compared  to  Zxx,3  =  495  —  j7969  using 
spectral  domain  integration.  For  T  =  0.1,  the  mutual  impedance  of  the  image  can 
be  calculated  by  numerically  double  integrating  Equation  (J.13a)  multiplied  by  the 
weighting  function,  resulting  in  ZIXi,  =  94.4  +  j'530.2,  compared  to  Zxx,t  =  95-f  j530 
using  spectral  domain  integration  by  subtracting  the  impedances  for  T  — >  00  and 
T  =  0.1.  Changing  the  parameters  to  W  —  0.03  and  kz  —  0.5 ku;  the  results 
are  Zxx<3  =  616.9  —  y‘35330  using  Equation  (1. 10)  and  Zxx,g  =  616  —  j33952  using 
spectral  domain  integration.  The  mutual  impedance  of  the  image  for  T  =  0.1  was 
ZXXi,  =  —89.2  +  j’560.3  numerically  and  Zxx>,  =  —89  +  j’560  by  spectral  domain 
integration. 

The  integrals  of  the  small  argument  approximations  of  the  Hankel  functions  used 
in  this  Appendix  are  summarized  in  Appendix  L. 
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Appendix  J 

Two-Dimensional  Traveling  Wave 
Electric  Currents 

This  appendix  presents  the  fields  of  a  two-dimensional  traveling  wave  electric  current 
given  by 

J(r')  =  J(f!)  e-f-'  (J.  1) 

in  a  two-dimensional  region,  R.  Thf  vector  potential  of  an  arbitrary  current  in  a 
volume  V  with  a  wave  number  k  is  calculated  using 

A(r)  =  h  / ,(r,)  e~r dv'  (J2) 

r 

where  r  —  |r  —  r'|.  Since  the  region  is  independent  of  z,  the  z  integration  of  Equa¬ 
tion  (J.2)  can  be  done  using 

dz'  =  -j  e-]k‘z  H[J\klP)  {J. 3) 

from  Appendix  K,  where 


r  =  y/p2  +  {z-z')2 

(J.4a) 

p  =  \/(*  —  x')2  +  (y  —  y')2 

( J  .4b) 

h  =  Jk2  -  kj. 

(J.4c) 

151 


The  resulting  equation  for  the  vector  potential  of  an  arbitrary  traveling  wave  electric 
current  in  a  ^-independent  region  is  given  by 

A  (p)  =  -J-  e I  Hu\ktp)  J(p')  dR'.  (J.  5) 

The  z-independent  component  of  the  current  can  be  written  in  rectangular  co¬ 
ordinates  as 

J(p')  =  xJx(p')  +  yjy{p')  +  Zjz{p')  (J.6) 

where  each  component  produces  a  vector  potential  of  the  same  polarization,  i.e., 


*-4(p') 

-»  *^x(p) 

(J.7a) 

yjyip') 

-♦  yK(p) 

( J.7b) 

zJz{p) 

-*  *A4(p). 

(J.7c) 

The  electric  and  magnetic  fields  are  calculated  from  the  individual  vector  potentials 
using  H  =  VxA  and  E  =  -juipA  +  V(V  •  A). 

The  fields  of  a  x  polarized  traveling  wave  current  are 


XT  A  -  dA*  - 

H  =  -jkzAxy  -  —  2 


E  =  —  Axfc2  + 


d2Ax  .  dAx  . 
*  +  a^-y v  - lk‘-^  z 


(J.8a) 

(J.8b) 


where  the  explicit  dependent  on  p  has  been  dropped  for  clarity.  Similarly,  the  fields 
of  a  y  polarized  traveling  wave  current  are 

- w  .  .  .  d Ay  „  /  ,  „  \ 


H  =  jkzAyx  + 


B  -  -| ^Ti+{A>kl 

use  oxuy  \ 


and  the  fields  of  a  z  polarized  traveling  wave  current  are 

dAz  .  dAz  . 

H  - 


E  =  —  I ~jkz^~  x  -  jk Z^1  y\  +  (fc2  -  k2)  Azz 


dr  J  *  8y 


(J.9a) 

(J.9b) 

(J.lOa) 

(J.lOb) 
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The  partial  derivatives  in  Equations  (J.8)-(J.10)  reduce  to  the  derivatives  of 


H„\k,p)  given  below: 

=  -k,  H?\k,p) 

ox  p 


(J.lla) 


%-Hl?\k,p)  =  -ktV—^-Hf\k,p) 
dy  p 


(J.llb) 


=  - - - —  -  H\2\ktp)  -  kt  H^\k'p ) 


(J.llc) 


-  l  if  -  x,y  -  (y  -  yy  uw 


-K\klP)  =  k, 


HY\klP) 


k^X  P~  H^(k,p)  (J.lld) 
P 


=  k(  H^k,p) 

or  ps 


-  k^-^l  H«\k,p).  (J.lle) 


For  example,  the  x  and  z  components  of  the  electric  field  produced  by  the 
traveling  wave  electric  current 


J(p')  =  x  JT(p)  e 


(J.12) 


in  a  region  R  are  given  by 


U  V 


2  if _ trWt 


{x  -  x'  f  -  (y  -  y'f  (2) 


HY\k(P)\  Jx{p)  dR'  (J.13a) 


E,(Jr)  =  -^7  /  ~~  H*\kP)  Mp)  dR'. 


(J.13b) 


The  x  and  z  components  of  the  electric  field  produced  by  the  traveling  wave  electric 


current 


J(p')  =  z  J2{p)  e 


(J.14) 
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in  a  region  R  are  given  by 


EX(JZ)  = 

Ez(Jz)  = 


—  —  e  jkz2  j -  H[2\kp)  Jx(p')  dR'  (J.15a) 

n  p 

e~jkzZJ  ^2)(M  Jz{p')  dR'.  (J.15b) 

H 
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Appendix  K 

Spherical  Wave  Spectrum  of  a 
Cylindrical  Traveling  Wave 


This  appendix  develops  a  spherical  wave  spectrum  of  a  cylindrical  traveling  wave 
starting  with  the  cylindrical  wave  spectrum  of  a  spherical  wave  given  by 

a—j  Ar-v/ p2 +{z— z')3  1  _ 

.  =  —  /  //,' 2 1  ( p vV  -  u,1 (K.l) 

+  U  -  2kL 

which  can  be  developed  using  Equations  (5-136),  (5-137),  and  (5-138)  of  [115,  p. 
244]  and  is  given  in  abbreviated  form  in  Equation  (5-139)  and  in  Equation  (6.614.4) 
of  [141].  Let 


/(*)  = 


g -Jky/p2+(z-z')’ 

\/p2  +  (z  -  Z')2 


and 


F(u>)  =  %H^\py/k2  -  u>2)e~JUZ 
3 


Now,  using  the  Fourier  transforms 

i  00 

/(*)  = 

—  oc 
oc 

F(u)  =  J  f(z)e~}“!dz 

—  OC 

the  final  relationship  is  obtained: 

-  e-*"  =  /  .  -  - 

j  L  \[? + (*  - 


,-JV* 


dz. 


(K- 2) 

(K.  3) 

(K.4a) 

(K.4b) 

-5) 
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Equation  (10)  of  [139]  incorporates  a  similar  result  with  7  =  —  jui. 
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Appendix  L 

Small  Argument  Approximations 
for  Integrals  of  Hankel  Functions 


This  appendix  presents  small  argument  approximations  for  several  integrals  of  Han¬ 


kel  functions  multiplied  by  other  functions,  using 

H^\z)  *  1  —  j—  In  ^rz  (L.la) 

IT  l 

H\%)  *  I  +  (L.lb) 

Z  TTZ 


where  7  =  1.781. 


—  a 


x'()  dx'  %  2a 


(a-*)ln  y|* 


a) 


~  k 

+  (a  4-  x)  In  —  |x  -f  a|  —  2a 


(L-2) 


/ 

—  a 


H[2\k\x  —  x'|)  dx'  %  fea  +  j 


irk 


1 

x  -f  a 


(1.3) 


j  j  Hf,2,(k\x  -  x'\)  dx'  tu(i)  dx  =s  2o  -  (-^  +  jln2  +  Imkaj  (LA) 


II 

—a  —a 


1 


x  —  X' 


H\2)(k\x  - 


x'|)  dx'  iu(x)  dx  2s  ka  —  j 


8 

irka 


In  2 


(L.b) 
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where 


^V(a  —  i)  for  0  <  x  <  a 
^f(a  +  x)  for  —  a  <  x  <  0. 


(L- 6) 
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Appendix  M 

Limitations  of  the  Reciprocity 
Theorem  in  Two-Dimensions 


This  appendix  examines  the  limitations  of  the  reciprocity  theorem  in  two-dimensions 
and  shows  a  practical  example  where  the  reciprocity  theorem  is  not  valid  [142].  First, 
we  will  derive  the  reciprocity  theorem  following  the  classical  approach  [115,  sec.  3-8]. 

Consider  two  time-harmonic  sources  a  and  b  with  electric  and  magnetic  current 
densities  ( J“,  Ma)  and  ( J6,  M6),  respectively,  radiating  in  a  region  containing  a  linear 
isotropic  inhomogeneous  medium  (/i,  e).  Source  o  or  b  radiating  alone  in  the  region 
produces  the  fields  (Ea,Ha)  or  (Eb,H6),  where  the  fields  are  related  by 

V  x  H°  =  jweEa  4-  J°  (M.la) 

-  V  X  Ea  =  ju>/zHa  +  M°  (M.lb) 

or 

V  X  H6  =  ju>eEb  +  3b  (M.2a) 

-V  X  Efe  =  ju>tiHb  +  Mb  (M.2b) 

respectively. 

Equations  (M.l)  and  (M.2)  are  combined  using  the  vector  identity 

V*(AxB)  =  B*V  xA  —  A*V  xB  (M.3) 
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resulting  in 


V  •  (Eb  X  H“)  =  H°  •  V  X  E6  -  Eb  •  V  x  H° 


=  -jup H°  •  Hb  -  H°  •  Mb  -  ju>eEb  •  Ea  -  Eb  •  Ja  (M.4) 


V  •  (Ea  X  Hb)  =  Hb  •  V  x  Ea  -  Ea  •  V  x  Hb 

=  —ju)fiU.b  •  Ha  -  Hb  •  Ma  -  ju>eEa  •  Eb  -  Ea  •  Jb  (M.5) 


Subtracting  Equation  (M.5)  from  (M.4)  yields 

V  •  (Eb  x  Ha  -  Ea  X  Hb)  =  Hb  •  Ma  -  Ha  •  Mb  +  Ea  •  Jb  -  Eb  •  Ja  (M.6) 

which  is  the  point  form  of  the  Lorentz  reciprocity  theorem. 

Applying  the  divergence  theorem  to  Equation  (M.6)  yields 

j H  (Eb  X  H“  -  Ea  X  Hb)  •  dS  = 

s 

jjj  (Hb  •  Ma  -  Ha  •  Mb  +  Ea  •  Jb  -  Eb  •  Ja)  dv  (M.7) 
v 

where  S  is  a  closed  surface  over  the  volume  V  and  dS  is  the  outward  normal  at  S. 
This  is  the  general  form  of  the  Lorentz  reciprocity  theorem  [131,  pg.  49]. 

The  sources  a  and  b  are  confined  to  volumes  Va  and  V&,  enclosed  by  surfaces  Sa 
and  5b,  respectively,  as  shown  in  Figure  M.l,  where  volume  Vo  is  the  source-free 
region,  enclosed  by  the  surface  E  +  Sa  +  Sb •  Since  (J°,Ma)  and  (Jb,Mb)  are  zero  in 
K,,  the  right-hand  side  of  Equation  (M.7)  is  zero  for  V  =  V0,  resulting  in  the.Lorentz 
reciprocity  theorem  for  a  source-free  region 


(Eb  X  H°  -  E°  X  Hb)  •  dS  =  0. 


(M.8) 


For  sources  of  finite  extent  in  all  three-dimensions,  the  surface  E  is  taken  to  be 
2.  rphere  cf  rndiu?  r  —*  oc.  Then  the  radiation  condition 


(Eb  X  H°  -  Ea  x  Hb)  •  dS  =  0 


(M.9) 
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is  used  to  reduce  Equation  (M.8)  to 


(E6  x  Ha  -  Ea  X  H6)  •  dS  =  -  <fj>  (E6  X  Ha  -  EQ  X  H6)  •  dS  (M.ltf) 

Sa  Sb 

which  is  often  written  as 

HI  (E6  •  J°  -  Hfc  •  Mj  dv  =  jjj  (Ea  •  Jb  -  Ha  •  M6)  dv  (M.ll) 

^  a  ^ 

using  Equation  (M.7). 

However,  for  sources  of  finite  radial  extent  in  two-dimensions,  the  source-free 
volume  Vo  is  a  cylinder  with  two  axial  voids  and  the  suface  £  is  the  outer  surface  of 
this  volume  and  can  be  decomposed  into  three  surfaces: 

£  =  Ep-.oo  4-  £2_oo  +  £2__oc  (M.12) 

The  radiation  condition  only  applies  for  the  surface  £p_oo,  so  the  reciprocity  theo¬ 
rems  of  Equations  (M.10)  and  (M.ll)  are  not  valid  for  general  sources  of  finite  radial 
extent. 

Next,  we  will  consider  under  which  conditions  the  surfaces  £*-.00  and  £z__0O  do 
not  contribute  to  Equation  (M.8).  This  contribution  is  given  by 

JJ  (E6  xH“-E“x  Hb)  'dS+JJ  (Eb  x  Ha  -  Ea  x  Hb)  •  dS.  (M.13) 

•  00  £*— *  — oc 

Th?  most  common  fields  in  two-dimensions  can  be  described  as  traveling  wave 
fields  where  the  only  z  variation  of  the  fields  is  The  two  limiting  cases  are 

TEz  and  TMz  fields,  where  kz  =  0,  and  TEM  fields,  where  kx  =  k  =  Since 

the  regions  of  interest  are  source-free,  i.e.,  the  surfaces  £*_oo  aQd  £r__O0,  the  fields 
can  be  represented  as  the  sum  of  the  fields  generated  by  the  magnetic  and  electric 
vector  potentials  [115,  sec.  3-12]  given  by 

A  =  z9{x,y)e~jk " 

F  =  z${x,y)e~jk «*. 
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(M.14a) 

(M.14b) 


The  resulting  fields  are  given  by 

V(V-A) 


E  —  —V  x  F  —  jui/iA  4- 


+- 


juje 


JU)€ 


.dv  .dV  . 

-jk,x- - jkzy~  +z 

ox  dy 


(k2  -  k2)  *J  |  e~*kzX  (M.15a) 


H  =  V  X  A  —  juie F  + 

«.d*  . 

XTy-VT*) 


V(V-F) 

juH 


+  ■ 


JvfJ- 


f)Ob  z  \  1 

a®  -  +  z  (fc2  -  fc2)  $  J  (M.15b) 

First,  individually  consider  the  surface  integrals  of  Equation  (M.13)  point-wise 

(E6  X  Ha  -  E°  X  H6)  •  z  (M.16) 

where  z  is  from  dS  =  zdxdy.  Using  Equation  (M.15),  the  first  term  can  be  expanded 
in  rectangular  coordinates  as 


(E‘  X  H“)  •  i  =  I 

[[* 

[  \  dy  u)fi  dx 


dy  ujc  dx 
d*a  k°  d$a 


j  V  j  e  Jk‘*  X 


\  dx  wc  dy 

)*-(?*  5? )>)■--  «•»> 


-{( 


) 


d$b  dya  d$b  dVa 
dy  dx  dx  dy 

klfd&cW  d$bd$a\ 
w/x  \  dy  dy  ^  dx  dx  ) 

:( 


^  kb  ld'Vbd'Sta  +  d^dy0' 


u;e  \  dx  dx 


dy  dy 


Wt  (W  w  _  w;\  1 

k2  \  dx  dy  dy  dx  )  j 


(M.18) 


163 


The  second  term  may  be  obtained  in  a  similar  fashion.  Then  Equation  (M.16) 
reduces  to 


(E6  x  Ha  -  E“  x  Hb)  z  =  | 

/  _  Kkz\  _  d$b  dVa 

\  k2  J  \  dy  dx  dx  dy 

8$a  d$bd$a\ 


<9£°\  (dVbd$a  d*bd$aY 
dy  )  \  dx  dy  dy  dx) 

\  i  (d*bdva  a^6^a\ll  /w 

J  a >e[dx  dx  +  dy  dy  jj}‘  (M'19) 


dx  dx 


Equation  (M.19)  is  zero  for  TEM  fields,  where  kz  =  k.  Using  Equation  (M.19), 
Equation  (M.13)  is  zero  for  traveling  wave  fields  where  k°  =  —kb,  and  as  a  spe¬ 
cial  case,  TEz  and  TM^  fields  where  k°  =  kb  =  0.  Under  limited  special  cases, 
Equation  (M.19)  may  be  zero  for  general  traveling  waves. 


Microstrip  Transmission  Line 

The  method  of  moments  solution  of  a  microstrip  transmission  line  is  one  case  where 
the  limitations  of  the  reciprocity  theorem  are  important.  This  problem  consists  of 
finding  the  electric  field  of  a  traveling  wave  current  when  calculating  the  method  of 
moments  impedance  matrix.  This  is  one  case  where  the  Galerkin  method  of  moments 
impedance  matrix  is  not  symmetric.  To  avoid  the  complexity  of  the  microstrip 
Green’s  function,  consider  the  mutual  impedance  between  the  currents  xJx(x)e~ikzZ 
and  zJ2(x)e~ik*z  in  free  space  using  Equations  (J.13b)  and  (J.15a) 
ir/2  [  .  M/2  f 

Zxz  =  -  /  -3-^~  j  ~^-H\2\k\x-x'\)Jz{x')dx'  Jx(x)dx  (M.20a) 

-M/2  -M/2 

MY*  .,  ,  U/2  f  _ 

z*r  =  -  /  “1^7  /  _  lgi2)(felg  ~  *'\)Mx)dz  Jz{x')dx  (M.20b) 

-M  /2  —  M/2 

where  kt  =  \Jk 2  -  k2.  In  this  special  case  Zxz  =  —  Zzx.  Clearly,  for  a  microstrip 
line  the  relationship  would  be  more  complicated,  but  definitely  Zxz  ^  Zzx,  except  in 
limited  special  cases. 
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